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In this Issue 

In 1965. V.W. Cooley and J.W. Tukey published their landmark "algorithm 
for the machine calculation of complex Fourier senes." which quickly became 
known as the fast Fourier transform, or FFT. The FFT reduced from hours 
to seconds the time required for a digital computer to calculate the frequency 
spectrum of a signal. Hewlett-Packard's first FFT analyzer, the HP 5450A 
Fourier Analyzer, made its debut in 1970. This general-purpose instrument 
could sample input signals between dc and 25 kHz, store the sample values 
in its memory, and later compute frequency spectra or perform other math- 
ematical manipulations of the data. 
Over the years, new twists and refinements have made the FFT ever faster, and advancing 
technology has increased processing speeds by a combination of faster digital circuits and parallel 
processing techniques such as pipelining. In this issue you'll find the design story of a new FFT 
analyzer that uses the latest technology, plus a few of its designers' own ideas, to sample signals 
between dc and 100 kHz and compute frequency spectra almost in real time — as many as ten 
spectra per second, fast enough that the user can watch the displayed spectrum change as 
changes occur in the system being analyzed. The HP 3561 A Dynamic Signal Analyzer is designed 
especially for electronic applications such as spectral analysis, network analysis, and waveform 
recording, for mechanical vibration applications such as predictive maintenance, dynamic balanc- 
ing, and spectral analysis, and for acoustic applications such as noise reduction and speech 
signal processing. Our cover illustrates one use of the HP 3561 As near-real-time processing, 
The plot shown is variously known as a spectral map, a stack plot, or a waterfall diagram. It's a 
series of successive spectra displayed simultaneously to show the results of a change of some 
parameter such as the speed of a motor. In this case the system is a fluid-film bearing (like the 
bearings on your car's crankshaft) and the spectral map shows the phenomenon called oil whirl 
becoming oil-whip instability as the shaft speed increases. For predictive machinery maintenance, 
the HP 3561 A is used just as a physician uses an electrocardiogram. A machine's vibration 
spectrum is compared with a normal one taken earlier to reveal changes that might signal an 
impending breakdown because of wear or internal damage. 

The HP 3561 A s capabilities are discussed in the article on page 4. On page 6 you'll find 
descriptions of some of its applications. Its hardware and software designs are described in the 
articles on pages 12 and 17, and the special aspects of its FFT implementation are outlined on 
page 20. The article on page 28 tells how the HP 3561 A s digital filters process so much data 
so quickly. 

Here's an invitation for HP end users and OEMs. Most of our articles are about the design of 
HP products and the contributions made by these products and designs. We'd like to publish an 
occasional article describing how HP hardware and/or software contributions have helped you 
make a contribution in your field. We're looking for original technical material that's interesting to 
more than just a few readers. If you have an idea for an article you'd like to write, send an outline 
to Editor, Hewlett-Packard Journal, 3000 Hanover Street, Palo Alto, California 94304 U.S.A. 



-fl.P. Dolan 



What's Ahead 

The design of a new line of general-purpose fiber optic test instruments will be featured in the 
January 1985 issue. The HP 8150A Optical Signal Source produces modulated, calibrated light 
signals for stimulating fiber optic devices. The HP 81 51 A Optical Pulse Power Meter measures 
various parameters of light signals, including peak levels, mesial power level, and extinction ratio. 
A third instrument, the HP 81 51 9A Optical Receiver, is a stand-alone optical-to-electrical transducer 
that gets users into fiber optic testing at the lowest possible cost. 
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Versatile Instrument Simplifies Dynamic 
Signal Analysis at Low Frequencies 

Analysis of low-frequency signals has many uses in 
electronic design, vibration studies, and acoustic 
measurements. This easy-to-use analyzer covers the range 
from 125 fiHz to 100 kHz and displays the data in several 
useful formats. 



by James S. Epstein 

REAL-TIME ANALYSIS of low-frequency electrical, 
mechanical, or acoustic phenomena can be of great 
value to design and maintenance personnel. Until 
the advent of digital signal processing, such analysis was 
difficult if not impossible to do. As the complexity and 
performance of digital ICs increased, it became possible to 
study low-frequency phenomena in more detail. 

Using fast Fourier transform (FFT, see box on page 20) 
digital analysis techniques, HP introduced several low-fre- 
quency signal analyzers in the 1970s— the HP 5450A 
Fourier Analyzer, the HP 5420A Digital Signal Analyzer, 
and the HP 3582 A Spectrum Analyzer. Combining many 
of the features of these earlier instruments with new capa- 
bilities made possible by VLSI circuit technology, the HP 
3561A Dynamic Signal Analyzer (Fig. I) is the latest 
member of HP's family of FFT analyzers. 
The HP 3561 A is designed to provide high-performance 



dynamic measurement and analysis capability in a small 
portable package. With the increased computational capa- 
bility that exists in today's microprocessors, many features 
can be now incorporated into an instrument that previously 
required a companion controller or computer. Analysis of 
dynamic signals has only recently been available because 
it awaited the development of suitable high-speed digital 
signal processing capability. This capability is incorpo- 
rated in the HP 3561A in the form of multiple micropro- 
cessors and custom digital integrated circuits. 

Conventional swept-frequency analyzers cannot analyze 
a dynamic signal very well because the signal may be 
changing during the sweep time. In dynamic signal analyz- 
ers, the signal is first digitized quickly and then analyzed 
after the data is captured. The analysis of the data may 
consist of examining the time waveform directly or examin- 
ing the frequency spectrum at evenly spaced intervals or 



Fig. 1. The HP 3561 A Dynamic 
Signal Analyzer combines high- 
performance hardware and ver- 
satile menu-driven software in a 
portable package for measure- 
ment and analysis ol low-lre- 
quency waveforms in electronic, 
acoustic, and vibration applica- 
tions Besides spectral analysis, 
the HP 3561 A can perform net- 
work measurements and time do- 
main waveform recording. 
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in Vi- or full-octave intervals. In addition, mathematical 
operations can be performed on the data to look at differ- 
ences, ratios, and other computed results. After the data 
has been captured and the desired results computed, the 
next step is presenting the data for display. The HP 3561A 
has several display formats, some of which allow up to 
sixty spectra to be viewed simultaneously. 

Key Features 

The HP 3561A can measure signals over the frequency 
range from 125 fiHz to 100 kHz. This covers the range of 
slow mechanical and geophysical vibrations to well above 
the 20- kHz audio band. One of the important specifications 
of a dynamic signal analyzer is its ability to resolve signals 
that are closely spaced in frequency. This allows a user to 
detect close-in sidebands on audio signals or vibrations 
caused by two different sources that are close in frequency. 
The HP 3561 A can analyze signals with a resolution of f>40 
/j.Hz over its entire 100-kHz range. 

Another important consideration in dynamic signal 
analysis is measuring small signals in the presence of large 
signals, such as the distortion of a sine input. A measure 
of an analyzer's ability to do this is called dynamic range; 
it specifies the maximum difference in signal levels that 
can be measured and guarantees that any internally gener- 
ated distortion and spurious products will be below the 
minimum measurable level. In the HP 3561A the dynamic 
range is 80 dB. which means that a primary signal and a 
secondary signal one ten-thousandth of the amplitude of 
the primary signal can be measured simultaneously. The 
accuracy of the measurement is also important. Using an 
internal calibration process, the HP 3561A has an 
amplitude accuracy of ±0.15 dB over the entire HP class 
B environmental range (0 to 55°C, <95% relative humidity 
at 40°C, and altitudes <4570 m). 

Applications 

The applications for a dynamic signal analyzer are very 
diverse. However, the HP 3561 A Dynamic Signal Analyzer 
is specially designed to cover three areas of application. 

The first area is general electronic measurements in the 
audio frequency range. One common measurement is the 
distortion of a pure tone. The HP 3561A can measure dis- 
tortion products thai are 80 dB below the fundamental and 
has a built-in marker function that computes the total dis- 
tortion of up to 20 harmonics when the marker is placed 
on the fundamental. Another electronic measurement is 
frequency response. The HP 3561A has an internal source 
that, combined with the internal trace arithmetic, can he 
used to measure both amplitude and phase as functions of 
frequency. General spectral analysis is another measure- 
ment done in the audio band to determine the exact fre- 
quency and amplitude of an input signal or to determine 
the frequency and amplitude of a spurious input. Transient 
analysis is another measurement often needed in the elec- 
tronic and mechanical areas. The 40,960-word buffer in 
the HP 3561A allows 0.8 second worth of data to be cap- 
tured for signals in the 0-to-20-kHz frequency range. 

Another major application area for dynamic signal 
analysis is measuring mechanical vibration using a trans- 
ducer that converts either displacement, velocity, or accel- 



eration into an electrical signal that can be measured by 
the analyzer. The analysis of vibration can be extremely 
useful in the maintenance of machinery. By examining the 
vibration in the frequency domain, small vibrations caused 
by worn-out parts can be detected even in the presence of 
large vibrations that are normal to the operation. In this 
manner, the machine can be shut down for repair only 
when necessary, saving costly down time. In addition to 
detecting when a problem has occurred, potential causes 
can be identified. (See box on page 6 for an example). 

A third application area for which the HP 3561A is par- 
ticularly well suited is acoustic measurements. Measuring 
acoustic signals requires special techniques. The human 
ear has its own frequency response, which must be con- 
sidered when making measurements. An internal A- 
weighting filter allows the analyzer to duplicate the weight- 
ing applied by the human ear. thereby allowing the noise 
levels measured to be compared to their levels as perceived 
by a human observer. Another important consideration in 
acoustic analysis is frequency resolution. The human ear 
cannot detect two frequencies separated by a fixed differ- 
ence at high frequencies as well as it can at low frequencies. 
Its response is logarithmic and is best characterized by 
using V:i-octave analysis bands that become wider as the 
center frequency increases. These analysis bands have been 
used for a long time to measure acoustic signals and many 
noise specifications are written in terms of Va-octave inter- 
vals as defined by ANSI specifications. The HP 3561A can 
display its results in 'A-octave, full-octave, or linear fre- 
quency increments to allow the user to determine which 
mode is most useful for a particular application. 

Friendliness 

With all of the capability that can be packed into an 
instrument such as the HP 3561 A, it could be difficult to 
use these features if the instrument were not specifically 
designed for ease of use. A measurement should be easy 
to make without requiring a lot of setup to get accurate and 
quick results. This is made possible in the HP 3561 A by 
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Dynamic Signal Analysis 
for Machinery Maintenance 



Using the HP 3561 A Dynamic Signal Analyzer to analyze the 
vibrational signature of machinery can allow identification of prob- 
lem spots before they cause machinery shutdown. In the time 
domain, vibrations from many sources in a machine can add to 
create complicated waveforms When viewed in the frequency 
domain, however, individual problems such as imbalances, gear- 
mesh anomalies, or bearing defects can usually be identified 
(Fig. 1) By saving the spectra of a healthy machine in the HP 



a damaged gear as shown by the excessive response at Ihe 
gear-mesh frequencies (Fig. 2b). 

In addition, the HP 3561 A has a spectral map display format 
that is particularly useful in servicing rotating machinery These 
spectral maps can be created by adding a third dimension to 
spectral measurements, such as time, r/min, or load The HP 
3561 A's map mode can be used to identify changes or trends 
in the spectra of machines, circuits, or other devices In the exam- 
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Fig. 2. (a) Baseline spectrum for a machine, (b) Later spectrum for machine of (a) showing 
the excessive response at the gear-mesh frequencies caused by a damaged gear. 



3561 A's optional one-megabit bubble memory, a set of baseline 
measurements is formed that can be used for later comparison. 
Fig 2a shows an example of a normal baseline measurement 
for a machine The same machine at a later time has developed 



pie shown in Fig 5 on page 10. a collection of measurements 
made during a machine runup forms a spectral map, illustrating 
changes in machine vibration with an increase in rotational 
speed 
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several features coupled with softkey menus. An autorange 
algorithm selects the correct range to place the input signal 
within 4 dB of full scale. Autocalibration is periodically 
performed so that the measurement accuracy is maintained 
within the = 0.15-dB specification. The first measurement 
covers the entire 100-kHz frequency range so that all fre- 
quencies can be examined. Then the user can zoom in on 
a particular range of frequencies by using the HP 356lA's 
marker function to identify the center frequency and press- 
ing the DEFINE CENTER softkey in the menu obtained by 
pressing the FREQuency key on the front panel. After the 
center frequency is selected, the desired analysis band can 
be selected in the FREQ softkey menu by pressing the DEFINE 
span softkey. 

Several special marker functions are provided to make 
the analysis of the measurement easier. The relative marker 
function allows any signal to be selected as the reference 
for either amplitude or frequency or both. When this marker 
is used, all measurements are scaled relative to the selected 
reference. This function is particularly useful when 
measuring how far the amplitude of a spurious signal is 
below thai of the desired signal. 

Another special marker function will mark up to twenty 
harmonics on the display screen with a special grid after 
the fundamental frequency is entered by the user (Fig. 2). 
In addition to marking the harmonic components on the 
screen, the total harmonic distortion is computed for the 
harmonics marked. A similar marker function is provided 
for marking sidebands and computing Ihe relalive power 
in Ihe sidebands with respect to the carrier. 

A very common type of analysis is computing the power 
in a specified frequency band. This is done when analyzing 
noise power in different frequency regions. The band- 
power marker function allows the user to position right 
and left markers to denote the band of interest and have 



the power in the band computed. 

Saving Instrument States 

Some specialized measurements may require many of 
the HP 3561 As special features to be selected. When this 
is the case, the storage and recall of instrument states in 
nonvolatile memory can be a very handy feature. Once the 
instrument is set up as desired, its state can be stored in 
one of six locations. If the optional bubble memory is in- 
stalled, up to 127 instrument states can be saved. At a later 
time, say at a remote measurement site, the desired state 
can be quickly recalled when it is time to make the measure- 
ment. This is not only convenient, but can help prevent 
errors in setting up the instrument. 

Narrowband Measurements 

The narrowband measurement is the basic analysis mode 
of Ihe HP 3561 A. This mode provides a frequency spectrum 
of 400 spectral lines equally spaced over the chosen 
analysis band. The analysis band can be centered anywhere 
in the instrument's frequency range from dc to 100 kHz 
with 0.25-Hz resolution and the frequency span for a 
zoomed measurement can be as small as 256 mHz. This 
gives a resolution of 640 fiHz anywhere in the spectrum 
up to 100 kHz. In addition, when making a baseband mea- 
surement (not zooming around a center frequency), the 
analysis band, or frequency span, can be reduced to 10.24 
mHz. which gives a maximum resolution of 25.6 nHz. 
Speed. Speed was an important consideration in designing 
the hardware to make narrowband measurements. When 
making adjustments to electronic equipment, such as re- 
ducing distortion components in a device under lest, Ihe 
display must be updated fast enough to give useful feedback 
to the person doing the adjustments. When averaging many 
spectra, the total measurement time is greatly dependent 
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upon the time required to generate each spectrum. The HP 
3561A can update a spectral display 10 times per second 
and can average up to 20 spectra per second. See the art Icle 
on page 12 for design details about how this performance 
was achieved. 

Windows. The shape of the analysis band is determined 
by the window applied to the data. In dynamic signal 
analysis, the band shape is selected by windowing — mul- 
tiplying the data in the time domain by a weighting func- 
tion. This is equivalent to filtering the data in the frequency 
domain. Just as it is useful to have different passband 
shapes in a conventional swept-frequency analyzer, it is 
also useful to have different windows in a dynamic signal 
analyzer. The HP 3561A offers a choice of four window 
functions — flattop, uniform, Hanning, and exponential. 
The flattop window gives maximum amplitude accuracy 
at the expense of frequency resolution. Its bandwidth is 
approximately four times that of the uniform window, 
which weights all of the time samples equally- The uniform 



window is used for signals that are periodic or are fully 
contained within the analysis band (such as a transient or 
an output of one of the HP 3561 A's internal noise sources), 
and it provides maximum frequency resolution. 

The Hanning window provides a good compromise be- 
tween frequency resolution and amplitude accuracy. The 
window flatness is within 1.5 dB, compared to 0.01 dB for 
the flattop window. Its bandwidth, however, is only 1.5 
times the maximum available resolution. 

The exponential window is useful when measuring 
transients. Coupled with the HP 356lA - s internal impulse 
noise source, this window allows the HP 3561A to make 
impulse response measurements. The shape of the expo- 
nential window can be adjusted to fit the data that has 
been collected. 

Measuring Phase. Another distinction between dynamic 
signal analyzers and conventional swept-frequency analyz- 
ers is the ability of the former to measure phase. There are 
also applications for measuring phase with respect to a 
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trigger, such as when balancing rotating machinery. Here, 
the measured phase is used to determine where in the 
rotation cycle an out-of-balance condition occurs. By using 
the HP 356lA's internal calibrator, the measured phase 
value is corrected for an accuracy of ±2 degrees up to 10 
kHz and ±10 degrees from 10 kHz to 100 kHz. 
Averaging. Averaging can be used to increase the useful- 
ness of the narrowband measurement. Several types of av- 
eraging are available. When measuring the spectral charac- 
teristics of random data, the accuracy of the results is sub- 
ject to the statistical properties of the data. One method of 
reducing the variance of the results is to perform rms or 
power averaging of the data. Another type of averaging that 
is useful for reducing the amount of noise on the signal is 
time averaging. This type of averaging requires a synchroni- 
zation signal, which can be either an external trigger or a 
part of the signal being measured. Two other types of av- 
eraging provided by the HP 3561 A are rms exponential and 
peak hold. The first is a type of rms averaging that weights 
earlier averages less than the current measurement. This 
allows a continuous update of the display while averaging. 
The second is used when the maximum value in each fre- 
quency bin is desired for a series of measurements. 
Triggering. Triggering is important when measuring phase 
and in time averaging Therefore, the HP 3561 A has several 
triggering modes. The source of the trigger can be the input 
signal, an external TTL signal, the internal source, or the 
internal trigger (free run). In addition, the instrument can 
be armed for a single trigger or armed to respond each time 
a trigger is received. Pretriggering of up to 8 records in 
narrowband mode or 40 records in time capture mode (dis- 
cussed later) is provided. Post-triggering of up to 1023 rec- 
ords can be selected. Pre- and post-triggering are very useful 
when a time difference exists between a known event and 
the data that is desired. When testing loudspeakers, for 
example, a time delay exists between the time the impulse 
excitation is generated by the instrument and when the 
impulse response is detected at the speaker. 
Simultaneous Display. The HP 3561A has the ability to 
display a full 100-kHz baseband spectrum and a zoom spec- 
trum simultaneously (Fig. 3). This capability is provided 
by the use of custom digital filters, which are discussed in 
the article on page 28. One way that this display format 
can be used is to have the full 100-kHz spectrum on the 
upper half of the display and the zoom spectrum on the 
bottom. Moving the marker to the desired signal on the 
upper display and pressing the CTR FREQ softkey in the 
marker menu will move the center frequency of the zoom 
spectrum to that of the desired signal. This can be per- 
formed repeatedly if desired to zoom in on several different 
signals for closer examination. 

Time Capture Measurements 

Time waveform capture can be used for transient analysis 
and the analysis of nonstationary signals such as speech. 
It also can be used to increase the real-time rate of the 
measurement to 100 kHz until the buffer is filled. 

The major difference between narrowband and time cap- 
ture measurements is that in time capture mode the entire 
time record is captured and stored in memory before any 
processing is done. If the bubble memory option is in- 



stalled, the time record can be retained even while the 
instrument is turned off. The time record can be defined 
to be from 1 to 40 records of 1024 samples each. The full 
triggering capability of the HP 3561 A can be used with 
time capture mode so that up to 40 records of pretrigger 
data can be taken, allowing the examination of data before 
the occurrence of a particular event. The ability to preserve 
the original time data before processing allows the center 
frequency and span for zoom analysis to be selected after 
collecting the data. This can not only save a lot of data 
collection time, but may be the only way to collect data 
when the input signal is of a transient nature. 

In time capture mode, the data is first collected within 
the frequency range specified for collection. If the data is 
collected in baseband mode, it can be translated to any 
center frequency up to the maximum frequency collected 
and zoomed up to a factor equal to the number of records 
collected. 

The postcollection translation can occur as many times 
as desired. In addition, the type of window applied to the 
data can also be changed to allow trade-offs between fre- 
quency resolution and accuracy. The postcollection pro- 
cessing is done by selecting the time capture data memory 
as the input to the digital filters. Because the data has 
already been collected, it can be run through the digital 
filters at maximum speed. Different start points in the data 
can be selected so that only desired data is analyzed and 
changes in spectral parameters can be examined at different 
points in the record. In addition, an average magnitude 
display can be selected to examine the average of all the 
data in the time capture memory at one time. 

When the time capture mode is combined with other HP 
3561 A features, very useful measurements can be made. 
For example, the average magnitude display can be used 
in the upper half of the display to show the full frequency 
range of the captured data, and the lower half of the display 
can then be set up to scroll through the captured data with 
or without zoom. When a user desires to look at the entire 
time buffer at one time, the spectral map display mode 
(discussed later) can be used. Whether scrolling through 
the data or mapping it, an increment can be set up to specify 
the percent of the data record to be incremented for the 
next display. In this manner, a spectral map can be dis- 
played so that each new trace in the map is offset by one 
to 1024 data points. This can be important when examining 
the decay of selected frequencies. 

Third- and Full-Octave Analysis 

For certain applications, most notably in the acoustic 
field, proportional frequency resolution is desirable. In this 
type of spectrum analysis, the signal is analyzed by a bank 
of filters with logarithmically spaced center frequencies 
and bandvvidths proportional to the center frequency. The 
most widely used analysis of this type is at Vj-octave and 
full-octave intervals. The HP 3561 A computes the Vii-octave 
and full-octave spectra digitally by combining narrowband 
frequency data taken on three separate analysis spans. To 
do this, the custom digital fillers are set up to collect data 
at three sample rates: f B , f s /10, and f,/100. The three data 
buffers are then transformed using the FFT. and a weighted 
sum of adjacent FFT filter bins is computed to form each 
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Vi-octave or full-octave filler bin. 

The frequency response of the octave filter shapes is 
determined by the windowing function used, the number 
of FFT bins summed, and the weighting function applied 
to the FFT bins. In the HP 3561 A. the Hanning window 
function is used for all of the '/i-octave and full-octave 
filter shapes. The number of bins summed varies with the 
center frequency of the octave filter. The weights applied 
to the adjacent FFT bins equal unity except for the first 
and last bins. 

The octave filter shapes are designed to meet ANSI and 
IEC standards. Because of the number of different octave 
spans allowed in the HP 3561A, 168 special filter shapes 
had to be designed and verified. Since this would have 
been a formidable task if it had been done manually, a 
computer program was written to find the filter weights. 
This program optimized the filter weights based on the 
noise equivalent bandwidth and the geometric mean center 
frequency specifications, which were the severest design 
constraints. A closed-form solution was used for the noise 
equivalent bandwidth and the geometric mean center fre- 
quency was optimized by using a binary search. The pro- 
gram completed the design task after a month of unattended 
operation on an IIP 9000 Series 200 Computer! 

Front Panel 

The front panel on the HP 3561 A is specifically designed 
for ease of use and ruggedness. The user interface with the 
instrument is through keys grouped by function and as- 
sociated softkey menus that direct the user through the 



available options. The major functional areas are Measure- 
ment. Trigger, HP-IB, Display, Marker, Instrument State, 
and Input. These groups are clearly indicated on the front 
panel so that the user knows which keys to press for any 
desired function. In addition to physical grouping, the dif- 
ferent types of keys are color coded so that related functions 
are closely associated. Within each functional group are 
specific keys that either perform a specific action or bring 
up softkey menus for further selections (Fig. 4). The softkey 
menus displayed along the right edge of the CRT screen 
show only the choices that are available for the mode 
selected. 

The keyboard is made up of a new conductive rubber 
key set that is custom molded for the HP 356lA's front 
panel. This keyboard allows multiple colors and very sharp 
lettering to be molded into the keys. The keys have a tactile 
feedback that indicates that a key has been pressed. In 
addition to their attractive appearance, the conductive rub- 
ber keys provide a ruggedness that is commensurate with 
the portability of the instrument. These rubber keys are 
very difficult to damage, even if objects are dropped on 
them. Furthermore, they provide resistance to dust and 
moisture and are very easy to install, reducing the assembly 
time and the cost of the instrument. 



Softkeys 

The softkey menu interface allows the user to be guided 
through the selections required for a setup. When a key, 
such as the MODE key in the Measurement group, is pressed, 
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Fig. 5. Spectral map display. Up to 60 magnitude traces can be displayed simultaneously lor 
studies of changes in the spectrum as a lunction of time, r/min. etc. These maps can be saved 
and recalled later from the optional nonvolatile bubble memory 
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the associated choices then available to the user are dis- 
played along the right edge of the display [see Fig. 1). For 
the MODE key. the resulting choices are NARROW BAND. 
THIRD OCTAVE. FULL OCTAVE. TIME CAPTURE. EXT SAMP ON 
OFF, DEFINE PULS/REV. and TEST SELECT. This is a good 
representation of the different types of softkeys available. 
The first four softkeys make up a group of selections, only 
one of which is allowed at a time. When the selection is 
made, the associated softkey label is intensified so that it 
is obvious which softkey is presently in effect. The next 
type of softkey is the on/off softkey or two-state softkey. 
This softkey allows the user the choice of having the labeled 
function either on or off. For example, the EXT SAMP ON OFF 
softkey allows the user to supply an external sample clock 
instead of using the internal one. When an on/off softkey 
is pressed, it toggles between on and off with its current 
state indicated by either the ON or OFF in its label being 
intensified. Another type of softkey is typified by the DE- 
FINE PULS REV softkey. This softkey requires a numeric 
entry, as do all softkeys that have the word DEFINE in their 
labels. When the label for this type of softkey is intensified, 
any numeric entries made will be for the parameter indi- 
cated on its label. If more than one DEFINE softkey label is 
available on a menu, only one will be intensified at a time. 
The units for the entry will be displayed on the softkey 
menu after any numeric key of the front panel's keypad or 
the intensified DEFINE softkey is pressed. 

The final type of softkey available is the select softkey 
as represented by the TEST SELECT softkey. When a softkey 
with SELECT in its label is pressed, another menu is dis- 
played with more selections. (See article on page 17 for 
more detail about the HP 356lA's software.) 

Display Formats 

Once the data is collected, the next step is to display it 
in a manner that besl allows the user to analyze the results. 
Because of the many different types of measurements pos- 
sible. I he type of analysis performed will vary depending 
upon the application. For this reason, several different dis- 
play modes are provided. We have already discussed the 
upper/lower display format in which two separate traces 
of any type can be displayed simultaneously, one above 
the other. There are many uses for ihis and il is the format 
that the instrument powers up in. Some examples are look- 
ing at magnitude and phase at the same time, observing 
the time waveform along with its magnitude or phase spec- 
trum, or viewing the full 100-kHz baseband spectrum along 
with a zoom spectrum. In the time capture mode, viewing 
a compressed time record and a magnitude spectrum to- 
gether allows the user to see which part of the large lime 
record is currently being processed. 

A single display format is also provided so that a display 
can be expanded to lake up the entire screen when only 
one trace is necessary to display the measurement. Another 
display format that is very useful is the front/back display. 
This display also allows any two traces to be viewed simul- 
taneously. However, they are bolh expanded to fill the 
entire screen with one positioned behind the other. The 
back trace is half the brightness of the front trace so that 
the two traces can be easily compared. The front and back 
traces can be switched by pressing the NEXT TRACE key. 



The most versatile display format is the spectral map 
display (Fig. 5) obtained by pressing the MAG MAP CLEAR 
softkey in the display format menu. Up to 60 magnitude 
traces can be simultaneously displayed, one behind the 
other, with a hidden-line algorithm applied so that the dis- 
play is more viewable. The traces can be positioned either 
directly behind each other or slightly offset along the hori- 
zontal axis so that data that may be blocked by the hidden- 
line algorithm can be seen. Through the use of the NEXT 
TRACE key and the Marker group keys, the full accuracy of 
the instrument is available while viewing any one of the 
60 traces. The NEXT TRACE key allows stepping through the 
map with the selected trace intensified. The marker func- 
tion then reads out the x and y values for the location of 
the marker to full precision. If it is desirable to expand the 
selected trace to fill the screen, a toggle key allows the user 
to switch between the spectral map and the selected trace. 
The data in the map is preserved in the HP 3561 A's memory 
until a new map is generated by pressing the MAG MAP 
CLEAR softkey again. If the bubble memory option is in- 
stalled, the entire map can be stored in nonvolatile memory 
for recall later. 

Portability 

Dynamic signal analysis applications are not always con- 
fined to the lab bench. These measurements are often made 
in remote locations such as on aircraft or at remote com- 
munications sites. Vibration analysis can be performed on 
machines that are located on shop floors, in power plants, 
in chemical factories, or in submarines. Hence, the HP 
3561 A is packaged much the same as a portable oscillo- 
scope. It weighs 33 pounds. A convenient carrying pouch 
is attached to the top of the instrument to store cables, 
manuals, and other equipment. 

Portability implies more than just a package; if a lot of 
additional equipment is required, portability is decreased. 
For example, when the instrument is carried to a measure- 
ment site, it may be desirable to bring the data back to a 
central computer for further analysis or storage. The op- 
tional bubble memory for the HP 3561A allows up to 127 
traces to be stored and retained even when the power is 
turned off. The traces can be recalled at any time for further 
analysis or transfer to another computer. 

When making vibration measurements, some kind of 
transducer is necessary. To reduce the number of transduc- 
ers that have to be carried around with the instrument, 
inexpensive integrated circuit piezoeleclric (ICP) ac- 
celerometers can be left attached to the machine under test. 
The current source required to power them is built into 
the HP 3561 A as a standard feature. 
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Hardware Design for a Dynamic Signal 



Analyzer 

by James S. Epstein, Glenn R. Engel, Donald R. Hiller, Glen L. Purdy, Jr., Bryan C. Hoog, 
and Eric J. Wicklund 



DESIGNING THE HARDWARE for the HP 3561A 
Dynamic Signal Analyzer involved the constraints 
of portability, adequate shielding of the signal pro- 
cessing circuits against electrical noise, a compact power 
supply, high-speed digital circuitry, and an accurate clear 
display of data and calculated results. Supporting and 
simplifying some portions of this hardware design are an 
extensive softkey menu structure and a sophisticated 
operating system (see article on page 17). 

Hardware Design 

A block diagram of the HP 3561 A is shown in Fig. 1. 
The front end performs the signal conditioning required 
to scale and filter the data so that it can be digitized by the 
analog-to-digital converter (ADC). After the signal is dig- 
itized, it is further processed by the digital filters to provide 
frequency translation and increased resolution. The main 
processor controls the entire instrument to perform the 
tasks required and communicates with the display to show 
the results of the measurement. The power supply provides 
the power to run all the internal circuitry and the non- 
volatile memory stores traces, setups, and time records for 
further analysis. An optional bubble memory can be added 
to extend the available memory. 

Front End 

The input amplifier circuitry (Fig. 2) in the HP 3561A 



provides a floating, shielded input connection for rejection 
of common mode and ground loop induced signals. A 
grounding switch is included for chassis-referenced mea- 
surements. The signal path does not become ground refer- 
enced until it leaves the ADC and crosses an isolation trans- 
former. Therefore, all circuitry through the ADC runs on 
separate floating power supplies. 

The input signal passes through a 0-. 20-. or 40-dl3 at- 
tenuator (depending on range) into a JFET input amplifier 
stage. Four fixed-gain amplifiers, each preceded by a pro- 
grammable attenuator, provide a total gain of 38 to - 40 dB 
in 1-dB steps. This gives the HP 3561 A a range of full-scale 
input values from 27 to -51 dBVrms in 1-dBVrms steps. 

The full-scale rms output voltage level of the input 
amplifier is 228 mV. Gain is distributed through the 
amplifier chain to maintain a reasonable signal level for 
best noise figure and distortion performance. The resulting 
noise floor specification is -141 dBV/VHz for frequencies 
lessthan 1 kHz and - 150 dBV/VHz above 1 kHz. Distortion 
is less than —80 dB relative to full scale. Internal offset 
voltages are corrected by an eight-bit DAC-controlled auto- 
zero circuit, which sums an opposing offset voltage into the 
last amplifier stage. A binary search algorithm is used by 
the firmware to set the D-to-A converter. 

Following the input amplifier is a passive nine-pole anti- 
aliasing filter with a rolloff above 100 kHz. The signal is 
then buffered, sampled, and digitized by the ADC. 
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Fig. 1. Hardware block diagram ol the HP 3561 A. 
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Additional circuitry includes the 4-mA current source 
for use with external integrated circuit piezoelectric (ICP) 
accelerometers. An internal A-weighting filter can be selected 
to provide an ANSI standard simulated "human ear" transfer 
function in the amplifier path. 

A-to-D Converter 

Requirements of the ADC are set by the bandwidth and 
dynamic range specifications of the instrument. The anti- 
aliasing filter technology used allows band widths up to 
approximately 0.4 times the sample frequency, which is 
fairly close to the theoretical Nyquist limit of 0.5. This 
results in a sample rate of 256,000 samples per second for 
the instrument bandwidth of 100 kHz. Dynamic range con- 
siderations imply a resolution of 13 bits and maximum 
nonlinearity of 14 LSB (least-significant bit). 

The HP 3561 A"s ADC achieves its high accuracy by using 
a technique known as residue conversion. The final result 
is formed by performing successive passes or partial A-to-D 
conversions of increasing accuracy and decreasing range. 
Each new pass starts by converting the ADC output (the 
result of the previous passes) to an analog signal using a 
DAC and then subtracting it from the ADC input. The result- 
ing residue represents the errors in the previous approxima- 
tion. The residue is then converted with a resolution and 
accuracy much greater than for the previous pass. The re- 
sult is added to the previous approximation of the ADC 
input to produce a much more accurate approximation. 

A partial block diagram of the ADC is shown in Fig. 3. 
This two-pass converter uses the same local ADC (an ADC 
within an ADC) to perform the conversions for each pass. 
During the first pass the ADC input signal is routed to the 
local ADC through an 8:1 attenuator. The result of the local 
ADC conversion is the first-pass approximation of the input 
and is stored in latch LI. This first-pass approximation is 
the input to the upper eight bits of the 13-bit DAC, which 
converts with a linearity equivalent to 'A LSB of 13 bits. 

The second pass begins as the first-pass approximation 



is latched and the DAC begins to settle. The DAC output 
is subtracted from the ADC input signal and the result is 
the residue representing the errors in the first-pass approx- 
imation. The residue is multiplied by four by an amplifier 
stage before the local ADC performs the second -pass ap- 
proximation. The gain ratio between the first and second 
pass is 4/(1/8) or 32. so the second pass resolves with 32 
times greater accuracy, but has a range 32 times less. The 
second-pass approximation is then added to the first-pass 
approximation in adder Al, but is shifted down five bits 
to account for the gain ratio of 32 between the first pass 
and the second pass. This final result is stored in latch L2 
before being transmitted to the digital filters. 

The most important aspect of a residue converter is the 
cancellation of first-pass errors. If an error E is made in the 
first pass, the DAC is loaded with the correct approximation 
to the input plus E. This is subtracted from the input so 
that the residue is only -E. This is then converted and 
added to the first-pass approximation to produce the final 
result so that -E, the residue, cancels the E term contrib- 
uted by the first-pass approximation. 

An important parameter in a residue converter is the 
overlap between passes. The HP 3561A's ADC has three 
bits of overlap between the first and second passes, mean- 
ing that the upper three bits of the second-pass approxima- 
tion have the same binary weight as the lower three bits 
of the first-pass approximation. The overlap is a measure 
of the magnitude of residue that can be converted by the 
second pass without overloading. If the first-pass approx- 
imation were perfect, then only quantization errors would 
be present in the first-pass approximation and no overlap 
would be required. 

The significant sources of linearity errors in this design 
are the DAC errors, the ADC errors on the second pass, and 
the analog gain matching between the first and second pass 
(how exactly the gain ratio of 32 is held). The gain matching 
determines how much the first-pass errors contribute to 
the nonlinearity of the converter. Since, in the best case, 
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the first-pass approximation still has quantization errors 
(the conversion is only eight bits, not 13), the gain matching 
is always important- 
Most ADCs (and DACs) have particularly nasty forms of 
nonlinearity. Common forms are step and impulse voltage 
transfer function errors. A step error is a constant error that 
is only added to the output when the input is in a certain 
range. This can occur when the analog weight associated 
with a particular bit is slightly in error. A common occur- 
rence is when the MSB (most-significant bit) of a DAC is 
in error. In this case a constant error is present whenever 
the input is greater than half scale. An impulse error is an 
error that is very localized. These nonlinearities result in 
distortion products that do not decrease (and can even rise) 
as the input level is lowered. This type of behavior is not 
desirable in a spectrum analyzer. One technique particu- 
larly effective for reducing this behavior is dithering (also 
discussed in the article on page 28). 

Dithering is the addition of an extra signal to the input 
to "smear out" nonlinearities. This causes any particular 
input level to correspond to different points on the voltage 
transfer function at different limes. In the long term the 
net transfer function at a particular input level becomes 
the weighted average over the neighborhood of that level 
in the original transfer function. The weighting function 
is the probability distribution function for the dither signal. 
Another way of looking at the effect of dithering is that 
harmonic products tend to be converted into intermoduln- 
tion products between the input and dither signals when 
the nonlinearities are localized as typical ADC and DAC 
nonlinearities are. This intermodulation effect implies that 
the dither signal should be broadband enough to distribute 
these products into many frequency bins so that the abso- 
lute magnitude at any frequency is much lower than if the 
products were distributed over only a small group of fre- 
quencies. This is equivalent to randomizing the errors so 
that they are not correlated with the input signal. 

A natural choice for a broadband dither signal is random 
noise. However, dithering only smooths out nonlinearities 
over a range the size of the dither signal, so distortion 
improvement is negligible for input levels greater than the 
dither level. This implies that the dither signal must be of 
a greater amplitude than the input signal to receive any 



distortion improvement. One way to increase the level of 
dither while not increasing the noise floor is to subtract 
the dither signal after the conversion. However, this rules 
out many sources of dither signals since the dither signal 
must be easily available in digital form. Hence, the use of 
pseudorandom noise is a natural choice. 

The HP 3561A's ADC uses this form of dithering in a 
special way. For the dithering signal lo be available both 
as an analog signal and as a digital signal implies that some 
sort of additional data conversion must be present, usually 
a DAC. It is possible to exploit the residue converter topol- 
ogy and achieve the same result without requiring another 
DAC. Fig. 4 shows the ADC with the dithering circuitry 
added. 

Since the first-to-second-pass gain match is not perfect, 
it is necessary to dither the first-pass conversion to ran- 
domize first-pass nonlinearity errors. Here the dithering 
signal Nf looks just like any other first-pass error and so 
is subtracted out in the combination of the first-pass and 
second-pass conversions. Since this dither signal appears 
in the DAC input and hence the residue, the DAC and 
second-pass conversion are somewhat dithered. However, 
since the first pass quantizes Nf, the dither signal has only 
several values and is not able to smear out nonlinearities 
adequately. A more continuous range of dither values is 
provided by adding pseudorandom noise to the first-pass 
conversion using adder A2. This dither signal also acts like 
a first-pass error since it is added before latch Ll. The result 
is that all of the sources of nonlinearity are dithered without 
the use of an extra DAC or ADC. The only extra circuitry 
required is a simple analog noise source Nf (actually one 
bit of pseudorandom noise passed through an RC filter to 
produce a more uniform probability distribution), the 
pseudorandom noise source, and adder A2, which per- 
forms only carries since the noise source has a zero mean. 
The circuitry within the shaded area in Fig. 4 is realized 
with a gate array in the HP 3561A. keeping the added cost 
of this dithering scheme very low. 

Processor 

The main processor of the HP 3561A is an 8-MHz 68000 
microprocessor with 393,216 bytes of program ROM. This 
processor controls the hardware and the user interface and 
performs the display processing and much of the data pro- 
cessing (see Fig. 1). It interfaces to the other hardware sub- 
systems via two different 16-bit buses: the processor bus 
and the RAM bus. The processor bus interfaces the main 
processor to the standard CMOS and optional bubble non- 
volatile memories, the display hardware, the local oscil- 
lator (LO), the HP-IB (IEEE 488) interface, the source, and 
the keyboard. 

The HP 3561 A has 131 .072 bytes of dynamic RAM, which 
is accessed through the RAM bus. This RAM is organized 
as 64K 16-bit words, but its data can be accessed as bytes 
by the main processor. This RAM is used for the processor's 
stack and program variables and for all the data buffers 
needed during the data processing. The RAM bus is con- 
trolled by a bus arbiter and the dynamic RAM controller. It 
has DMA ports to the dedicated FFT processor and the 
outputs of the digital filter control chip. The digital filter 
control chip performs DMA accesses using two address 
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counters which control where in RAM the data is written. 
The RAM bus also has some control and status ports used 
by the main processor for interfacing to the FFT processor, 
the digital filter control chip, the digital filter DMA address 
counters, timing logic, and the interface to the HP 3561 A s 
front-end circuitry- and ADC. 

Nonvolatile Memory 

The standard nonvolatile memory in the HP 3561 A con- 
sists of 8K bytes of CMOS RAM. This RAM is used to store 
data for two traces and six instrument setups, with one of 
the setups being the last state before power-down. The 
CMOS RAM is backed up by a three-volt lithium battery, 
which allows it to retain data for up to five years. 

There is also a bubble memory option that can add a 
megabit of nonvolatile memory. This allows users to store 
extended time records, as well as more setups and traces. 
The bubble memory can be accessed by a simple processor 
interface and can transfer data at a rate of 12. 5K bytes/s. 
These speeds allow the user to retrieve an extended time 
trace of forty records in 0.7 second, worst case. 

Noise Source 

The noise source in the HP 3561A can generate three 
types of signals: periodic, impulse, and random. All three 
of the signals are band-limited and band-translated. Like 
the tracking oscillator in a swept analyzer, the noise source 
always tracks the selected frequency span. 

The impulse source has the added capability of being 
triggered from the HP 3561 A's front panel or via the HP-IB. 
and it can skip time records. These features are useful when 
transient inputs are required to settle before a new measure- 
ment is made. The impulse source can skip up to 32,766 
records before generating another impulse. 

Digital Filter 

The digital filter section is based on custom ICs de- 
veloped especially by HP for dynamic signal analysis. Con- 
trolled by the main processor, this section filters the dig- 
itized input waveform to form the appropriate data records 



for the FFT processor. See the article on page 28 for a 
detailed discussion of the design constraints and realiza- 
tion of this section. 

Display 

The HP 3561 A uses a raster-scanned magnetic-deflection 
CRT for display of measurement results. This is a departure 
from some other Hewlett-Packard signal analyzers, which 
use a vector-scanned electrostatic-deflection CRT display. 
The raster-scan technology offers significant advantages for 
the HP 3561 A. Probably the most important is compact 
size. With magnetic-deflection CRTs, large electron beam 
deflection angles are easily achieved. Because of this, mag- 
netic-deflection CRTs can be much shorter than electrosta- 
tic-deflection CRTs for a given viewing area. The space 
saved by using a magnetic-deflection display helped fit the 
HP 356lA's extensive measurement capability into an os- 
cilloscope-sized package. Another advantage of the raster- 
scan technology is low cost. The HP 3561 A display, includ- 
ing drive electronics, costs about as much as an electrostatic 
CRT by itself. 

Display quality was not sacrificed to gain these advan- 
tages. Because of their shorter electron beam path length, 
magnetic-deflection CRTs give bright, crisply focused dis- 
plays. Vector scanning is limited in the number of traces 
that can be displayed before flickering becomes annoying. 
Raster scanning does not have this problem. This allows 
the HP 3561 A to display measurement results in three-di- 
mensional form with up to 60 traces of 400 points each 
(see Fig. 5 on page 10). That's the equivalent of over 24.1)00 
vectors. Also, raster scanning consistently produces good- 
quality alphanumerics. Well annotated displays are an im- 
portant part of the friendliness of the HP 3561A. 

Display quality is further enhanced by scanning the raster 
vertically at a rate of 36,205 Hz with the vertical lines 
formed from left to right. Since most rapid trace variations 
occur vertically, this gives a smoother trace by reducing 
the jagged line appearance compared to a horizontally 
scanned raster. 

The digital portion of the display contains a memory 
control circuit, a special sync circuit, and a processor inter- 
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face. The display is bit mapped with a vertical resolution 
of 256 pixels and a horizontal resolution of 512 pixels. The 
pixel bit map is stored in two sets of eight 16K x l dynamic 
RAMs. One set is for full-intensity pixels and the other is 
for half-intensity pixels. Different intensities are used to 
give emphasis to specific: areas on the CRT. To allow the 
display to be synchronized to the power line frequency, a 
special sync circuit is used which does not require phase- 
locked loops or analog circuitry. The circuit free-runs dur- 
ing display updates and synchronizes during retrace inter- 
vals when minor timing variations are not noticeable. The 
addition of the synchronization circuit requires only one 
more flip-flop. 

The processor interface allows quick generation of dis- 
play graphics. It includes alphanumeric generation capabil- 
ity and simple line-drawing hardware that can draw verti- 
cal, horizontal. 45°. and dotted lines at 1.400.000 pixels 
per second. It also supports several block operations such 
as set. clear, and complement at 1 1 ,300,000 pixels per second. 
These hardware features significantly reduce the amount of 
software overhead required for a bit-mapped display. 

Power Supply 

Requirements for the HP 356lA*s power supply were 
that it be lightweight and compact for instrument portabil- 
ity and that it provide eight different outputs, four of which 
are isolated floating supplies. A switching supply design 
seemed to be the natural choice for the job. The potential 
disadvantage of using this type of supply is that it can be 
a source of EMI (electromagnetic interference). For this 
reason, an oul-of-band switching frequency of 128 kHz syn- 
chronized to the system clock is used. 

A half-bridge, off-line circuit topology was chosen, using 
500-volt VMOS switching transistors (Fig. 5). Control loop 
circuitry regulates the main 5V logic supply. The output 



voltage is normally under control of a voltage feedback 
loop, but during overcurrent conditions a separate loop 
takes control to provide foldback current limiting. Both 
loops have their own loop shaping networks. All other 
supplies are derived from other windings on the same trans- 
former and thus are slave preregulated by the regulation 
of the 5V supply. Each of these other supplies is then further 
regulated by a linear three-terminal regulator. This pro- 
vides additional ripple rejection and supply isolation and 
the preregulation maximizes the efficiency of the linear 
regulators. 

The power supply is protected by four separate protec- 
tion circuits. The four fault-detection circuits are floating, 
nonfloating. and 5V overvoltage comparators and a pri- 
mary-overcurrent-sense circuit. All protection circuit out- 
puts are monitored by a latching shutdown circuit that turns 
off the gate drive to the power VMOS transistors and indi- 
cates on one of four LEDs which fault condition occurred. 

Packaging and Shielding 

The HP 3561 A package is similar to that used for portable 
oscilloscopes. The circuit boards sit in a vertical card cage 
and plug into a motherboard. Since the display and the 
switching power supply were identified as potential EMI 
sources, they are both placed behind a shield that runs the 
length of the instrument. 

Suppression of residual spurious signals by more than 
-80 dB relative to full scale required shielding of both the 
display and input amplifier. Electromagnetic field frequen- 
cies of concern are the display's vertical scan rate of 36,205 
Hz and its harmonics. Electric fields caused by the high 
accelerating potentials for the CRT are attenuated by a 
grounded screen that covers the CRT face. High-impedance 
circuits in the HP 356lA's input amplifier are covered by 
a grounded aluminum shield can. 
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Varying magnetic fields are generated by the display's 
magnetic elements (flyback transformer and yoke) and by 
other ac high-current carrying conductors. These fields can 
generate currents in the instrument frame, causing local 
field radiation near the input amplifier. The initial field 
levels near the input amplifier were small (milligauss 
levels), but the performance specifications required addi- 
tional attenuation of more than 40 dB. Shield materials 
have lower permeabilities at these low field levels, making 
the task more difficult. To simplify the problem, field radi- 
ation is reduced at the flyback transformer by a shorted 
turn of aluminum wrapped around the transformer. Cir- 
culating currents in this loop create opposing fields (Lenz's 
law), to help cancel some of the radiation, which in turn 
reduces field levels at the input amplifier. Further attenu- 
ation is achieved by enclosing the input amplifier board 
in a MuMetal '" shield. The CRT assembly (including yoke) 
is enclosed in another MuMetal shield, both to reduce its 
radiation and to shield the CRT from external fields that 
can cause display distortion. Fields from the power sup- 
ply's magnetic components have frequencies at 128 kHz 
and its harmonics and thus are above the HP 3561 A's mea- 
surement baseband of 100 kHz (one of the primary reasons 
for switching the power supply at 128 kHz). Aliased spuri- 

MuMetai is a iraaemar* cl Aliegneny-Luaium Siee! Corporation 



ous supply radiation problems are eliminated by syn- 
chronizing the supply to the system clock, placing the 
aliased products at dc. Careful attention to grounding tech- 
niques also played an important part in the total solution. 
The resulting spurious signal levels are less than — 80 dB 
relative to full scale on the HP 3561 A's most sensitive range 
(-51 dBVnns). 
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Instrument Software for Dynamic Signal 
Analysis 

by Glenn R. Engel and Donald R. Hiller 



THE SOFTWARE for a microprocessor-controlled in- 
strument determines the instrument's usefulness 
and personality — how it interacts with the user and 
how easy it is for a new user to understand and use effec- 
tively. Because there are many variables in dynamic signal 
analysis measurements, the software for the HP 3561 A 
Dynamic Signal Analyzer is particularly important. Some 
of the key features provided by the HP 3561A's operating 
system are: 

■ Friendly user interface 

■ Autocalibration 

■ Overlapped processing for quick display updates 

■ Spectral map display 

■ Direct output to HP-IB (IEEE 488) printers and plotters 

■ Nonvolatile memory storage 

■ Self-tests and service tests. 

Architecture 

At the heart of the HP 356lA's software is a 4K-byte 
multitasking operating system that provides process syn- 



chronization and resource control. There are five operating 
system tasks in the HP 3561 A (see Fig. 1). 

■ I/O control 

■ Command processing 

■ Measurement/display 

■ Nonvolatile memory control 

■ Message handling. 

The I/O control process interprets commands from the 
HP-IB and front panel. This process checks syntax and 
parameter limits and does immediate execution of simple 
commands. An example of a simple command would be 
setting the marker position. Commands that cannot be exe- 
cuted immediately such as buffer dumps and nonvolatile 
memory operations are passed on to the command process 
for completion. The command process then coordinates 
the operations necessary to complete the command. This 
includes locking necessary resources, communicating with 
other processes to handle parts of the command, and then 
unlocking resources after command completion. 

The measurement/display process is responsible lor all 
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signal processing functions and display updates. Signal 
processing functions include data collection, fast Fourier 
transforms of data (see box on page 20). amplitude correc- 
tions, phase corrections, averaging, and data scaling. The 
measurement/display process takes advantage of over- 
lapped processing to increase data throughput. This pro- 
cess is shown in Fig. 2 and is discussed later. 

The nonvolatile memory process is responsible for all 
nonvolatile memory functions such as saving and recalling 
traces and instrument states. With the optional bubble 
memory installed, the nonvolatile memory process be- 
comes very important, because of the slower access time 
of bubble memory devices. By handling memory operations 
with a separate process, the HP 3561 A can be making useful 
measurements while traces or states are being saved or 
recalled from the bubble memory. 

The message process handles error messages and user 
warnings. This function is a separate process because mes- 
sages need to be prioritized and most warnings have a 
display timeout. By using the message process, other pro- 
cesses do not have to worry about what happens to a mes- 
sage once the message process knows about it. 

The operating system provides several functions to ma- 
nipulate semaphores. The primary request and release func- 
tions provide a synchronization mechanism for resource 
allocation and are used to control counting semaphores for 
message communication between processes. When a pro- 
cess requires a resource, it calls the operating system's 
request function. If the resource is available, the process 
continues execution. If not, the operating system suspends 
the process and does not allow it to run until the resource 
is available. A resource becomes available when the process 
using it calls the operating system release function. When 
a resource is released, the operating system checks for any 
processes waiting for the released resource and schedules 
them for execution. 

Friendly User Interface 

The software for the HP 3561 A's user interface was writ- 
ten to provide the user with as much feedback as possible 
concerning error conditions and instrument operation. For 
example, if the user enters an illegal value for a numeric 
input, the HP 3561A, instead of just beeping, responds by 
telling the user exactly what the minimum and maximum 
values can be. In addition, the software helps the user by 
giving suggested operation hints. For example, when the 
pseudorandom noise source is selected, the HP 3561A re- 
minds the user that the uniform window is the most likely 
window selection for that noise source. If the uniform win- 
dow is already chosen, no message is given. The friendly 
user interface extends to the HP-IB interface as well as 
front-panel operation. A diagnostic mode allows for easy 
debugging of HP-IB error conditions. This mode displays 
each character received from the HP-IB interface on the 
instrument's display and marks the exact location of syntax 
violations, which can save programmers considerable de- 
bugging time. 

Autocalibration 

Calibration in the HP 3561 A is a three-step process. First, 
the software characterizes the amplitude and phase re- 



sponse of the instrument's analog front end over the entire 
operating range. Second, the software computes the front- 
end correction for the frequency range currently selected 
for measurements. This range, or span, can be a small por- 
tion of the HP 356lA's baseband of 100 kHz. Third, the 
software simulates the digital filters to determine the effects 
of the filtering operation on the magnitude and phase of 
the input signal. This results in one correction that repre- 
sents any errors introduced by the HP 3561A between its 
input connector and the final digital representation of the 
input signal in RAM. 

The key to providing high amplitude accuracy for the 
HP 3561 A is directly related to its ability to measure its 
own magnitude and phase response properly. An accurate 
digital-state-to-analog-level (1-bit D-to-A) converter deter- 
mines the accuracy of the internal calibration signal; it is 
driven by a pseudorandom noise generator. The spectrum 
of this calibration signal consists of a number of frequencies 
(determined by the length of the pseudorandom sequence) 
thai are of nearly constant magnitude and whose magnitude 
and phase are theoretically exactly known. This allows 
rapid calibration cycles, since both magnitude and phase 
of all desired calibration frequencies can be determined 
with a single operation. 

One problem of this autocalibration scheme is slew-rate 
limiting of the analog version of the signal, since the ideal 
signal has step changes between the two levels. If the slew- 
rate limiting is symmetrical, the effect is that of passing 
the signal through a low-pass filter with an impulse re- 
sponse that has the width of the slew time. For high slew 
rates (which are present in the HP 356lA's calibration cir- 
cuit), this filtering causes negligible errors. On the other 
hand, if the slew rates are not symmetrical, the effect is 
that of adding a small signal that is the derivative of the 
pseudorandom sequence. This additional signal causes 
small but significant errors in the amplitudes of the 
pseudorandom noise frequencies and the resulting calibra- 
tion constants. Fortunately, these errors are nearly identical 
and of the same phase when an amplitude-inverted 
pseudorandom sequence is passed through the same cir- 
cuitry. This allows the error to be canceled by subtracting 
a measurement of the inverted sequence from a measure- 
ment of the noninverted sequence. Another way of looking 
at this situation is that the errors are even-order and can 
be cancelled by the same technique used in a differential 
transistor pair or a push-pull amplifier to reduce even-order 
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distortion. Differences in propagation delays in the digital 
circuitry driving the analog conversion circuitry also cause 
the same kind of errors and so are also canceled. 

By using pseudorandom noise autocalihralion, the HP 
3561 A is able to correct for antialiasing filter rolloff well 
past the cutoff frequency of 100 kHz. However, because of 
possible aliasing problems and decreased dynamic range, 
the accuracy of the HP 3561 A is not specified above 100 
kHz. The resulting accuracy specification is ±0.15 dB from 
0 to 100 kHz. 

Overlapped Processing 

Overlapped processing is possible in the HP 3561A 
thanks to a multiport RAM. This RAM is accessible from 
four different sources: the main processor, the slave FFT 
processor, and two DMA (direct memory access) channels. 
This allows the HP 3561 A to almost triple its real-time rate 
by doing three operations simultaneously. Fig. 2 illustrates 
the possibilities. While the main processor is manipulating 
record N for display, the FFT processor is processing record 
N + 1, and the DMA sequence is collecting record N + 2. 

Spectral Map Display 

The software for the HP 356lA's spectral map display 
(see Fig. 5 on page 10) is quite simple. This display format 
provides a three-dimensional representation of the analyzed 
input signal versus time. As new spectra are placed into 
the map, they must be drawn as if they were behind the 
previous spectra. This operation is similar to viewing a 
distant range of mountains from an airplane. Mountain 
ranges that are behind others are only visible if they are 
higher than the ranges in front of them. That is, higher 
mountains can obscure lower mountains behind them. The 
hidden-line algorithm used by the HP 3561A uses exactly 
this idea. When a new spectrum is put on the display it is 
compared point by point with the previous spectrum to 
decide whether it is "visible." This process is simplified 



for a raster scan display because the two spectra only need 
to be compared at a limited number of points — once every 
horizontal pixel position. No interpolation is needed be- 
tween points and only the largest value for each horizontal 
pixel position must be remembered for later spectra, 

HP-IB Plots and Printer Dumps 

HP-IB plots and printer dumps are supported on the HP 
3561A to allow the user to produce a hard-copy represen- 
tation of the display. Printer dumps are achieved by simply 
reading the bit-mapped display memory and sending this 
information to a printer that supports raster-type graphics. 

Plots are obtained by sending HP-GL plotter commands 
to a pen plotter. A high-quality copy of the HP 3561A 
display is obtained in this manner. Plotting the spectral 
map is more complex than generating the raster-display 
spectral map. When plotting the spectral map, the HP 
3561 A must keep a copy in its memory of the largest values 
plotted to aid in hidden-line removal. These values are 
kept with enough resolution to provide correct hidden-line 
removal down to the width of the plotting pen. HP-IB plots 
also allow the user to annotate specific marker values to 
show test results or merely to provide extra documentation. 
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FFT Implementation 



One of Ihe benefits of using an analyzer based on a fast Fourier 
transform (FFT) is measurement speed. The HP 3561A Dynamic 
Signal Analyzer uses the FFT algorithm to analyze signals with 
the equivalent of a parallel bank of 400 narrowband filters This 
simultaneous spectral analysis technique is inherently much fas- 
ter than using a swept-frequency spectrum analyzer that must 
cover the same frequency span serially, one frequency al a time 
Although not yet practical at RF and microwave frequencies be- 
cause of the fast FFT computation rate required, current technol- 
ogy does allow its use at low frequencies This speed improve- 
ment can be put to good use in many of the HP 3561 A's appli- 
cations The typical display update rate of five to ten spectra per 
second makes circuit adjustments easier A noisy signal can be 
averaged 1000 times in less than a minute in the fast averaging 
mode. 

These fast measurements require a fast FFT computation. In 
the HP 3561 A. the FFT computation is performed by a TMS320 
microprocessor. This processor has a reduced instruction set 
designed for digital signal processing applications. The proces- 
sor and its support chips occupy only 30 square inches of printed 
circuit board area. It computes a 512-point complex FFT in 35 
ms and a 1024-point real FFT in 42 ms, including the window 
multiplications. 

Speed 

Fast execution times are obtained by tailoring the FFT algorithm 
to the internal architecture of the TMS320. A modified radix-8, 
Sande-Tukey twiddle-factor algorithm' is used. With a slight re- 
arrangement of the multiplications, this algorithm fits nicely into 
the TMS320's sum-of-products architecture. A radix-8 algorithm 
was chosen for several reasons With this radix, bit-reversmg is 
replaced by digit-reversing, that is, groups of bits are reversed 
at one time, which is easier and only needs to be done for every 
eighth data address. Also, a radix-8 algorithm requires fewer of 
the slower I/O operations to the instrument's main RAM To reduce 
the effect of slow RAM accesses further, the algorithm performs 
arithmetic operations using the on-chip cache RAM while waiting 
for the I/O operations 

Noise and Accuracy 

Noise performance of the FFT algorithm is equally important 
for the HP 3561 A. A numerical algorithm is only as accurate as 
the number of bits used in the calculations Inaccuracies in cal- 
culating the FFT can degrade absolute amplitude accuracy and 
linearity and generate phantom spectral components These are 
all lumped into the category of "noise " In the HP 3561 A, the FFT 
noise performance can be no better than the RAM data bus 
width, which is 16 bits. Given a perfect algorithm (no errors), the 
maximum dynamic range would be approximately 96 dB (6 dB 
per bit) This is more than adequate to ensure the instrument 
dynamic range specification of 80 dB, including all scaling over- 
heads and specification margins 

The challenge is keeping a 16-bit FFT algorithm accurate to 
1 6 bits at us output The algorithm reads data words from main 
RAM. performs arithmetic operations, rounds the results and 
stores them back in main RAM This occurs once for each FFT 
level, and there are four or five levels per transform If a 16-bit 
data pomt is multiplied by a 16-bit coefficient, the result is 31 
bits long If this result is then summed with eight others just like 
it, three more bits are added and the result could overflow the 
TMS320's 32-bit accumulator Therefore, prescaling is necessary 

The HP 3561 A s FFT algorithm uses a simplified block scaling 
technique It is implemented by saving the magnitude of the 



largest data point observed for an FFT level The next level 
chooses to divide by one, two, four, or eight based on Ihe previous 
level's maximum, This scaling method gams much of the perfor- 
mance improvement of a true block floating-point algorithm with- 
out the added complexity or speed reductions. With this type of 
scaling, the instrument noise floor drops with the level of the 
input signal. 

Another important consideration is the nature of Ihe noise floor. 
A noise floor that is uniform versus frequency, or white, is pre- 
ferred. This means that the average power of the error in each 
of the output frequency bins is the same. Most of the bins are 
the result of successive multiply-accumulate operations, where 
many bits are discarded in the rounding operation. It is fairly 
easy to keep the error in these bins consistent. The problem bins 
are the ones that are the result of an accumulate only In these 
bins, very lew bits are thrown away during rounding. The error 
has only a few possible values, and there is an expected error, 
or bias, in the result. For example, if two 16-bit values are added 
and the 1 7-bit result is rounded back to 1 6 bits, the error is either 
zero or Vi LSB (least-significant bit). There is a bias of V* LSB. 
When these biased results are summed in following levels, the 
bias error adds up in a linear fashion. 

Generally, FFT noise floors have spur-like bumps located on 
the output bins with large expected errors A user might be 
tempted to interpret these as spurious components in Ihe input 
signal even though they are outside of the specified dynamic 
range To solve this problem, two types of rounding offsets are 
used One results in a positive bias, and the other results in a 
negative bias. By appropriately choosing between these when 
rounding, the bias errors cancel. 

The FFT limits amplitude accuracy in yet another way: the 
windowing operation. The input data block is multiplied by a 
tapering function, or window, that forces Ihe start and the end 
of the record to zero This reduces the effecls of leakage, which 
occurs when the signal is not periodic in the measurement inter- 
val. Windowing is often thought of in the frequency domain as 
the frequency response of the analysis filters. To make accurate 
amplitude measurements, the effective filter frequency response 
should be very flat over a frequency span of plus or minus half 
of a bin spacing The HP 3561 A's flattop window function is flat 
to ±0.003 dB, much flatter than those used in other FFT-based 
analyzers. This flatter window is needed to complement the im- 
proved amplitude measurement accuracy of the HP 3561 A. 

Circuitry 

The TMS320 hardware looks like a slave processor to the HP 
3561 A's 68000 main system processor The 68000 sets up a 
small table in main RAM that tells the TMS320 what type of FFT 
to perform, scaling restrictions, data buffer start addresses, and 
other setup information Then the TMS320 S reset line is set true 
and FFT computation begins. When finished. Ihe TMS320 writes 
the scale factors in main RAM and interrupts the 68000 Main 
RAM is shared among several processors and DMA ports. A bus 
arbitrator establishes priority and grants RAM accesses. 

Acknowledgments 

Ron Potter designed the improved flattop window. 

Reference 

1 EO Bngham, The Fast Fouriet Transloim Prenlice-Hail. 1974 

Bryan C. Hoog 
R&D Engineer 
Lake Stevens Instrument Division 



20 HEWLETT-PACKARD JOURNAL DECEMBER 1984 



© Copr. 1949-1998 Hewlett-Packard Co. 



HEWLETT-PACKARD JOURNAL 

Volume 35 January 1984 through December 1984 

Hewlett-Packard Company 3000 Hanover Street Palo Alto. California 94304 USA 
Hewlett-Packard Cenirai Mailing Dept . Van Heuven Goedharttaan 12V 1 1 81 KK Amstetveen The Netherlands 
Hewlett-Packard (Canada) Ltd 6877 Goreway Drive Mississauga Ontaro L4V 1 M8 Canada 
Yokogawa-Hewtert-Packard Lid Suginami-ku Tokyo 168 Japan 



ncex 



PART 1: Chronological Index 



January 1984 

Two High-Capacity Disc Drives. Kent Wilken 

A Command Language for Improved Disc Protocol , Douglas L. V'oigt 

Second-Generation Disc ReadAVrite Electronics. Robert M. Batey 

and James D. Becker 
Disc Drive Error Detection and Correction Using VLSI. Peter M. 

Galen 

Head Positioning in a Large Disc Drive. R. Frank Bell. Eric W. John- 
son. R. Keith Whitaker, and Roger V. Wilcox 
Mechanical Design of a Large Disc Drive, lames H. Smith 
High-Capacity Disc Drive Servomechanism Design. Stephen A. 
Edivards 

Speech Output for HP Series 80 Personal Computers, Loren M. 

KoehJer and Timothy C. Mackey 
Speech Output for HP 1000 and HP 3000 Computer Systems. 

Elizabeth R. Hue/tie and Jeffrey R. Murphy 

February 1984 

A New Series of High-Performance Real-Time Computers, Marlu 
E. Allan. Nancy Schoendorf, Craig B. Challerton. and Don M. 
Cross 

An Adaptable 1-MIPS Real-Time Computer, David A. Fotland, 

Lee S. Moncton, and Leslie E. Ne.ft 
Designing a Low-Cost 3-MrPS Computer, Donald A. Williamson, 

Steven C. Steps, and Bruce A. Thompson 
Floating-Point Chip Set Speeds Real-Time Computer Operation. 

William H. McAllister find John R. Carlson 
Comprehensive. Friendly Diagnostics Aid A-Series Troubleshoot- 
ing, Michael T. Winters und /ohn F. Shellon 
New Real-Time Executive Supports Large Programs and Multiple 

Users. Douglas O, Hartman, Steven R. Kusmer, Elizabeth A. 

Clark, Douglas V. Larson, and Billy Chu 
New Software Increases Capabilities of Logic Timing Analyzer, 

David L. Neuder 

March 1984 

A New 32-Bit VLSI Computer Family: Part II— Software, Michael 
V. Helrick und Michael L. Kolesar 

HP-UX: Implementation of UNIX on the HP 9000 Series 500 Com- 
puter System, Scoll W. V. Wang and Jeff B. Lindberg 

An Interactive Run-Time Compiler for Enhanced BASIC Language 
Performance, David M. Landers. Timothy W. Tillson, Jack D. 
Cooley. and Richard R. Rupp 

A Local Area Network for the HP 9000 Series 500 Computers. John 
J. Balza, if. Michael Wenzel, and fames L. Willils 

Data Communications for a 32-Bit Computer Workstation, Vin- 
cent C. Jones 

A General-Purpose Operating System Kernel fora 32-Bit Computer 
System. Dennis D. Georg, Benjamin D. Osocky, and Stephen D. 
Scheid 

The Design of a General-Purpose Multiple-Processor System, 
Benjamin D. Osecky, Dennis D. Georg, and Robert /. Bury 

An I/O Subsystem for a 32-Bit Computer Operating System, Robert 
M. Lenk, Charles E. Meur, /r.. and Marcel E. Meier 

Viewpoints — Coping with Prior Invention. Donald L. Hammond 

April 19H4 

Low-Dispersion Liquid Chromatography. Robert /. /onker and 

Gerard P. Rozing 
Identification and Quantitation of PTH Amino Acids, Bernd GJotZ 



and Rainer Schuster 
Design of the HP 1090 Control System, Herbert Wiederoder. 

Roland Martin, and Juergen Ziegler 
A New Solvent Delivery' System. Wolfgang Ceiger and Heinrich 

Vollmer. 

Automatic Liquid Chromatograph Injection and Sampling. Wolf- 
gang Kretz and Hans-Georg Hart) 

Mobile Phase Preheater Ensures Precise Control of LC Column 
Temperature. Helge Schrenker 

A Low-Cost LC Filterphotometric Detection System, Axel Wiese, 
Bernhnrd Dehmer. Thomas Dorr, and Giinter Hoschele 

A High-Speed Spectrophotometry LC Detector, /oachim Leyrer, 
Gunter E. Nill. Detlev Hadbaivnik, Giinter Hoschele. and /oachim 
Dieckmann 

New Technologies in the HP 1090 Liquid Chromatograph. Alfred 
Maute 

May 1984 

Putting a 32-Bit Computer System in a Desktop Workstation. 

lack L. Burkman. Robert L. Brooks. Ronald P. Dean. Paul F. 

Febvre. and Michael K. Boiven 
Color Graphics Display for an Engineering Workstation, Daniel 

G. Schmidt 

BASIC Language Graphics Subsystem for a 32-Bit Workstation, 
Kenneth W. Lewis, Alan D. Ward, and Xuan Bui 

I/O Features of Model 520 BASIC. Gary D. Fritz and Michael L. 
Kolesar 

A Compact. Reliable Power Supply for an Advanced Desktop Com- 
puter. Jack L. Burkman, Howell R. Felsenthal. Thomas O. 
Meyer, and Warren C. Prat! 

Compact 32-Bit System Processing Units, Kevin W. Allen, Paul 
C. Christo/anelli, Robert E. Kuseski. Ronald D. Larson. Duvid 
Maitland, and Larry /. Thayer 

June 1984 

A Parametric Test System for Accurate Measurement of Wafer- 
Stage ICs. Yoh Narimatsu and Keiki Kanafuji 

Powerful Test System Software Provides Extensive Parametric 
Measurement Capability. Tukuo Bunno 

A High-Speed 1-MHz Capacitance/Conductance Meter for Measur- 
ing Semiconductor Parameters. Tomoyuki Akiyamu and Kenzo 
Ishiguro 

An Electronic Tool for Analyzing Software Performance. Gail E. 
Hamilton. Andrew /, Blasciak, Joseph A, Hawk, and Brett K. 
Carver 

Counter Module Simplifies Measurements on Complex Wave- 
forms. Donald /. Smith, Johnnie L. Hancock, and Thomas K. 
BohJey 

July 1984 

A New Handheld Computer for Technical Professionals, Suson I.. 
Wechsier 

Soft Configuration Enhances Flexibility of Handheld Computer 

Memory, Nathan Meyers 
Custom CMOS Architecture for a Handheld Computer. James P. 

Dickie 

Packaging the HP-71B Handheld Computer. Thomas B. Lindberg 
Moduli! Adds Curve-Fitting and Optimization Routines to the 1 171 IS. 

.Stanley M. BJdSCOM , Jr. and James A. Donnelly 
ROM Extends Numerical Function Set of Handheld Computer. 



(DECEMBER 1984 HEWLETT-PACKARD JOURNAL 21 



© Copr. 1949-1998 Hewlett-Packard Co. 



Laurence VV. Grodd and CharJes M. Pallon 
Plug-In Module Adds FORTH Language and Assembler lo a Hand- 
held Computer. Robert M. Miller 

August 1984 

Touchscreen Personal Computer Offers Ease of Use and Flexibility. 

Srinivas Sukumar 
Operating System and Firmware of the HP 150 Personal Computer, 

Laurie E. Pollero Wood and Charles H. Whelun 
The HP 150 Touchscreen: An Interactive User Input Device for a 

Personal Computer, Peter H. Stralon, Scoll R. McClelland, and 

Thomas L'. Kilbourn 
Applications Software for the Touchscreen Personal Computer, 

Peler S. Showman. Karl W. Pettis. Karlie /. Arkin. Jeffrey A. Spoel- 

stra. John Price. W. Bruce Culbertson. and Robert D. Shurlle/f. Jr. 
Hardware Design of the HP 150 Personal Computer. John E. Watkins. 

Patricia A. Brown, George Szeman, and Susan E. Carrie 
Software Graphics in the HP 150 

Personal Computer Printer Is User Installable, Joseph D. Barbera 
A Standard Keyboard Family for HP Computer Products, Lorenzo 
Dunn and Michael R. Perkins 

September 1984 

Transmission Impairment Measuring Set Simplifies Testing of 
Complex Voice and Data Circuits, David fl. Novolny, Jeffrey 
Tomberlin, Charles P. Hill. James P. Quan. Gordon A. Jensen, 
and Jerry D. Morris 

TIMS Mechanical Design 

Weight, Size, and Noise Impact Power Supply and Display Design 
Master/Slave TIMS Operation Improves Productivity, Teresa L.Reh 
How Master/Slave Mode Works 

Testing the TIMs. Allan W. Dodge, Scott S. Neal, and Kurt R. 
Goldsmith 

Semiconductor Research Corporation: A Perspective on Coopera- 
tive Research, Richard A. Lucie 
A Hyphenation Algorithm for HPWord, Paul H. Smil 
Designing Software for the International Market, Heather Wilson 
and Michael J. Shaw 

October 1984 

The HP 3065 Board Test Family: A System Overview. Thomas 

H. Fay and John E. McDermid 
HP Q-STAR 



Confirmation-Diagnostics 

Automatic Test Program Generation for Digital Board Testing, 

Robert E. Balliew 
Board Test Connection Terminology 

Digital Subsystem for a Board Test System. Matthew L. Snook 

and Michael A. Teska 
Digital Test Throughput 

Safeguarding Devices Against Stress Caused by In-Circuil Test- 
ing. Vance R. Harwood 

Extensive Library Simplifies Digital Board Test Setup, Handy W. 
ilolmberg 

An Interpreter-Based Board Test Programming Environment. 

Mark A. Malhieu 
Testing for Short-Circuit Failures. T, Michael Hendricks 
Reducing Errors in Automated Analog In-Circuil Test Genera- 
tion, John E. McDermid 

November 1984 

An Advanced 5-Hz-to-200-MHz Network Analyzer, Robert A. 

Wide and Jerry W. Daniels 
User-Defined Vector Math Expands Measurement Capabilities 
A Broadband Two-Port S-ParameterTest Sot. William M. Spaulding 
An ADC for a Network Analyzer Receiver. Alan J. Baker 
An Industrial Workstation Terminal for Harsh Environments, Jean 

Bounaix, Jean-Claude Dureau, and Jacques Firdmann 
How Do You Describe Terminal Ruggedness? 
High-Quality. Dot-Matrix Impact Printer Family, Murk J. DiVittorio 
Custom IC Controls Dot-Matrix Printers, Thomas B. Pritchard and 

David S. Lee. 

December 1984 

Versatile Instrument Simplifies Dynamic Signal Analysis at 

Low Frequencies, James S. Epstein 
Dynamic Signal Analysis for Machinery Maintenance 
Hardware Design for a Dynamic Signal Analyzer, James S. Epstein, 

Glenn R. Engel, Donald R. Hiller, Glen L. Purdy, Jr., Bryan C. 

Hoog, and Eric J. Wicklund 
Instrument Software for Dynamic Signal Analysis. Glenn R. Engel 

and Donald R. Hiller 
FFT Implementation 

Custom Digital Filters for Dynamic Signal Analysis. Charles R. 
Panek and Steven F. Kator 



PART 2: Subject Index 



Subject Month 
A 

Acoustic measurements Dec. 
Active input and output 

circuits, TIMS Sept. 

Actuator, disc drive |an. 

ADC. dynamic signal analyzer Dec. 

ADC. receiver, network analyzer Nov. 

ADC. 16-bit. low-cost Apr. 
Algorithm, board test connections Oct. 

Algorithm. Booth Feb. 

Algorithm. Fletcher-Powell July 

Algorithm, hyphenation Sept. 

Algorithm, short circuit testing Oct. 

Algorithm, touchscreen scanning Aug. 
Alignment, automatic. 

disc drive head )an. 

Amino acid analysis Apr. 

Amplifier, prescaling, ADC Nov. 

Analog testing, in-circuit Oct. 



Analyzer, dynamic signal Dec. 

Analyzer, logic timing Feb. 

Analyzer, network Nov. 

Analyzer, software performance June 

Architecture. 32-bit processing units May 

Architecture, handheld computer July 

Architecture, parametric test system June 

Array redimensioning July 

Assembly language ROM July 
Autocalibration. dynamic signal 

analyzer Dec. 

Autoinjector, LC Apr. 

Autosampler. LC Apr. 

Averaging, time interval June 

B 

Baluns, broad frequency range Nov. 

BASIC programmable computer July 

BASIC. 3-D graphics May 

BASIC. Board Test Oct. 



BASIC. HP 9000 Model 520 Mar. 

BASIC. I/O, HP 9000 Model 520 May 

BASIC, multiprogramming May 

Bed-of-nails test fixture Oct. 
Bell-standard transmission 

measurements Sept. 

BIOS Aug. 

Bit-slice data acquisition processor Apr. 

Board testing, in-circuit Oct. 

Bond-wire heating, in-circuit testing Oct. 

Booster pump Apr. 

Booth algorithm Fab. 

Bridge, reflectometer, configurations Nov. 

Bubblsjal Mar. 

Buffer, circular, I/O May 

c 

Cable compensation, capacitance 

measurements June 

Cable connections, sealed Nov. 
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Cache memory 


Feb. 


Development, parallel 


Mar. 


H 




Cache memory, L'O 


Mar. 


Di-bit integrator 


(an. 


Head alignment, automatic, disc drive 


Jan. 


CALC mode. HP-71B 


luly 


Diagnostic control system 


Feb. 


Head positioning system, disc drive 


Jan. 


Capacitance measurement 




Diagnostics. 32-bit SPUs 


May 


High-pressure pump 


Apr. 


subsystem 


June 


Digital filter 


Dec. 


High-speed columns. LC 


Apr. 


Card file software, personal 




Digital subsystem, board test 


Oct. 


Holographic diffraction grating 


Apr. 


computer 


Aug. 


Digital testing, in-clrcult 


Oct. 


Horner's method, polynomial roots 


July 


Card reader/recorder, handpulled 


July 


Diode array detector 


Apr 


HPLC 


Apr. 


Case, sealed 


Nov. 


Disc drives, 404M-byte 


Ian. 


HPWord hyphenation algorithm 


Sept. 


Ceramic disc, rotary valve 


Apr. 


Dispersion, holographic grating 


Apr. 


HP-IB transformer 


Apr. 


Clock recovery, disc drive 


(an. 


Dispersion. LC 


Apr. 


HP-IL. interface module 


July 


CMOS, handheld computer system 


luly 


Display, color graphics 


May 


HP Q-STAR, network, software 


Oct. 


CMOS testing, latchup prevention 


Oct. 


Display, dynamic signal analyzer 


Dec. 


HP-UX operating system 


Mar. 


Code and data separation 


Feb. 


Display, tilt mechanism 


May 


HP-UX. corporate strategy- 


Mar. 


Code segment mapping 


Feb. 


Dithering, dynamic signal analysis 


Dec. 


Hyphenation algorithm 


Sept. 


Code, intermediate 


Mar 


Dot-matrix printer family 


Nov. 






Codes, disc drive 


Ihii. 


Duty cycle measurements 


June 


1 




Coding, linear predictive 


(an. 






IEEE floating-point standard 


luly 


Coherence breaking by phase 




E 




In-circuit device testing 


Oct. 


modulation 


June 


EMI design, rope gasket 


May 


Industrial workstation terminal 


Nov. 


Columns, liquid chromatography 


Apr. 


Environmental IP rating 


Nov. 


Infrared array touchscreen 


Aug. 


Command Set 80 (for disc drives) 


Ian. 


Environments, industrial. 




Injector, LC 


Apr. 


Communications processor (COM) 


Apr. 


terminal for 


Nov. 


Ink jet printing, thermal 


Mar. 


Communications units 


Apr. 


Error correcting memory 


Feb. 


Integral function 


July 


Compiler, run-lime, interactive 


Mar. 


Error correction, capacitance 




Integrals, nesting 


July 


Complex data type 


luly 


measurements 


June 


Interfaces. I/O. TIMS 


Sept. 


Complex mathematics 


luly 


Error detection and correction, 




Intermodulation distortion, 




Computer. 32-bit 


Mar- 


disc drive 


Jan. 


measurement circuit 


Sept. 




May 


Errors, analog in-circuit testing 


Oct. 


Interrupt handling 


Mar. 


Computer, handheld 


luly 


Events measurements 


June 


IP rating 


Nov. 


Computer, personal. LC controller 


Apr. 


Exchanges 


Aug. 


IPG-II. automatic test program 




Computer, touchscreen, personal 


Aug. 






generator 


Oct. 


Computers, HP 1000 A-Series 


Feb. 


F 




ISO substitution tables 


Sept. 


Configuration, automatic 


June 


Fast Fourier transform 


July 


I-to-V converter 


In ne 


Confirmation/Diagnostics 


Oct. 


Fasteners, tolerance design 


Dec. 


I/O drivers 


Mar. 


Connection modes, capacitance 




May 


I/O subsystem, 32-bit computer 


Mar. 


measurements 


June 


FFT analysis, low-frequency 


Dec. 


I/O. BASIC, examples 


May 


Connector, custom power transistor 


Mav 


File management system 


luly 


I/O, TRANSFER statement 


May 


Contexts St 'N oueralinc svstem 


Mar. 


File manager, personal computer 


Aug. 


I/O unified 


May 


Control svslem LC 


Apr- 


File security 


luly 






Cooling, desktop workstation 


NT IV 


Filter, digital 


Dec. 


J 




Cooling, power supply 


Mav 


Filterphotometric detector 


Apr. 


Junction temperature, in-circuit 




Cooperative research 


Sept. 


Financial calculator software, 




testing 


Oct. 


Counter, gated universal 


(una 


personal computer 


Aug. 






Co-injection molding 


Sept. 


Fire code 


Ian. 


K 




Crystal, quarlz. measurements 


Nov. 


Fletcher-Powell method 


July 


Keep/toggle vectors 


Oct. 


Curvi* fittiuu 


lulv 


Floating-point chip set 


Feb. 


Key redefinition 


luly 


Cuslom IC nriuter control 


Nov. 


Flow cell 


Apr. 


Keyboard family, computer products 


Aug. 


1 ii' li.il If**; r 




FM code, Manchester 


Jan. 


Keyboard, detached design 


May 


f iiclnm IP c rtiuititl filttr 

i tit I.' 1' II 1 1 '1 llllCl 


Dec. 


FORTH language ROM 


luly 


Keyboard, sealed 


Nov. 


C-HIGH RESOLN function 


June 


FORTH/BASIC conversion 


July 






C-t mpasurpmenls 


)une 


Fractional-N synthesizer board 


Nov. 


L 




C-\' mpa^u rcmiTil s 


June 


Fractional-N synthesizer, modified 


Sept. 


Laguerre's method, polynomial roots 


July 




Frequency measurements 


|une 


LAN 9000. local area networking 


Mar. 


D 




Front panel, virtual 


June 


Level measurements, TIMS 


Sept. 


Daisy-chain configuration 


July 






Library, digital test 


Oct. 


Damping. LC 


Apr. 


G 




Library, speech 


Jan. 


Data acquisition processor (DAP] 


Apr. 


Gated universal counter 


June 


Library, waveform measurement 


June 


Data circuit testing. Bell-standard 


Sept. 


Generic research 


Sept. 


Lightweight frame, TIMS 


Sept. 


Data communications, workstation 


Mar. 


Graphics software, personal 




Limit cycles, digital filter 


Dec. 


Data storage and retrieval, logic 




computer 


Aug. 


Linguistic rules 


Sept. 


analysis 


Feb. 


Graphics, 3-D, BASIC 


May 


Liquid chromatograph 


Apr. 


Datacnm drivers 


Aug. 


Graphics, asynchronous input 


May 


LO board, network analyzer 


Nov. 


Debugger, software 


Mar. 


Graphics, input device tracking 


May 


Localization, software 


Sept. 


Degradation table 


Sopt. 


Graphics, text 


May 


Logic development, timing analyzer 


Feb. 


Detector, touchscreen 


Aug. 


Grating, holographic diffraction 


Apr. 


Longitudinal balance measurement. 




Detectors, liquid chromatography 


Apr. 


Group delay, measurement 


Nov. 


TIMS 


Sept. 
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Low-dispersion liquid 

chromatography Apr. 

Low-pressure compliance Apr. 

Low-pressure pump Apr. 

M 

Machinery maintenance Dec. 

Management, project Mar. 

Manchester code. KM Jan. 

Marked events, logic analysis Feb. 

Mask design, filter IC Dec, 

Master/slave mode. TIMS Sept. 

Matrix calculations July 

Measurement, dc subsystem June 

Mechanical design, 32-bit SPUs May 

Mechanical design, disc drives Jan, 

Mechanical design, power supply May 

Mechanical design, TIMS Sept. 

MemoMaker Aug. 

Memory, overlay structure Feb. 

Memory, soft-configuration July 

Memory, virtual Mar. 

Mentors, industrial Sept. 

Metering pumps Apr. 

MFM code Jan. 

Microbore columns, LC Apr. 

Micromachine. HP 1000 A-Series Feb. 

Microparaphraser Feb. 

Mixer-multiplier, digital filter Dec. 

Mixing chamber Apr. 

Mobile phase Apr. 

Modeling, foam core board May 
Modular design, desktop workstation May 

MS-DOS operating system Aug. 

Multiple-processor system design Mar. 
Multipoint connections, industrial 

terminal Nov. 

Multiprogramming, BASIC May 

Multiuser facility Feb. 

Mute character table Sept. 

N 

Narrowband measurements, 

low-frequency Dec. 

Network analyzer Nov. 

Networking, board test Oct. 

Networking, local area Mar. 

NLQ, near letter quality printing Nov. 

NMOS-III ICs May 
Noise reduction, TIMS display and 

power supply Sept. 



Octave analysis, acoustic 

measurements Dec. 

Offset errors, digital filter Dec. 

Offset frequency board Nov. 

Opens test Oct. 
Operating system, dynamic signal 

analyzer Dec. 

Operating system. HP-UX Mar. 

Operating system. MS-DOS Aug. 

Operating system, RTE-A Feb. 

Operating system, SUN Mar. 

Optimization, example July 
Oscillator. 100-MHz lime base. 



phase modulated 
Oscilloscope measurement system 
Oven, crystal reference 
Oven, LC 

Overlapped processing 



Packaging, handheld computer 
Parametric testing, semiconductor 
Peak-to-average ratio measurements 
Period measurements 
Personal Applications Manager 
Personal computer, touchscreen 
Phantom short circuits, testing 
Phase modulation, random oscillator 
Phase shift measurements 
Phonetics, speech synthesis 
Pipeline shift registers. 

digital filter 
Pipelined data path 
Point-to-point connections, 

industrial terminal 
Polynomial root finder 
Postprocessing, logic analysis data 
Power supply design, dual 
Power supply, desktop computer 
Power supply, dynamic signal 

analyzer 
Precision machining 
Preheater, LC mobile phase 
Printer, thermal, personal computer 
Printers, dot matrix 
Printhead. carriage control 
Printwires, drive control 
Prior invention, coping with 
Probing pattern generator 
Processing units. 32-bit system 
Production tests. TIMS 
Programming, board test 
Propagation delay measurements 
PTH amino acid analysis 
Pulse width measurements 
Pumps, metering 
Pump, booster 
Purity, LC peak 



RACE program 
Rating, IP 

Read/write electronics, disc drive 
Real-time computers, 1-3 MIPS 
Receiver, digital, IF 
Receiver. TIMS 
Reference board, frequency 
Regulator, primary switching 
Reliability assessment tests. TIMS 
Removable disc pack, 404M-byte 
Research, semiconductor, generic. 

cooperative 
Rise/fall time measurements 
Rotary valve. LC 
RTE-A operating system 
Rugged industrial terminal 
Run-time compilation 



Safeguards, in-circuit test 



June 
June 

Nov. 
Apr. 
Dec. 



July 
June 
Sept. 
June 
Aug. 
Aug. 

Oct. 
June 
June 

Jan. 

Dec. 
Feb. 

Nov. 
July 
Feb. 

Sept. 
May 

Dec. 
Apr. 
Apr. 
Aug. 
Nov. 
Nov. 
Nov. 
Mar. 
June 
May 
Sept. 
Oct. 
|une 
Apr. 
June 
Apr. 
Apr. 
Apr. 



Sept. 
Nov. 

Jan. 

Feb. 
Nov. 
Sept. 
Nov. 

May 
Sept. 

Jan. 

Sept. 

fwon 

Apr. 
Feb. 
Nov. 
Mar. 



Oct. 



Sampler. LC Apr. 

Sampling, software measurements |une 

S-parameter measurements Nov. 

Sapphire metering pump pistons Apr. 

Scanned infrared array Aug. 

Scanner card, digital Oct. 

Sealing, industrial terminal Nov. 

Sector sparing Jan. 
Semiconductor parameter 

testing June 
Semiconductor Research 

Corporation Sept. 

Separation, code and data Feb. 

Shorts test Oct. 
Signal analysis, dynamic, 

low-frequency Dec. 

Signed digit encoding Feb. 

Softkey menus, network analyzer Nov 

Software graphics Aug 

Software localization (translation) Sept 

Software translation Mar. 

Software. C/C-V Meter June 

Software, digital analyzer Dec. 

Software, multiple-processor system Mar. 

Software, parametric test system June 

Software, performance analysis June 

Software, preserving investment Mar 

Software, speech output Jan. 

Software, touchscreen computer Aug. 

Software, waveform measurement June 

SOLVE function July 

Solvent delivery system. LC Apr. 

Spectra, LC Apr. 

Spectral map display Dec. 

Spectrophotometric detector Apr. 

Speech output, HP 1000 Computer Jan. 

Speech output, HP 3000 Computer Jan 
Speech output. Series 80 Computers Jan 
Spreadsheet program, personal 

computer Aug. 

Spring, micro Apr. 

Stationary phase Apr. 

Statistics, software analyses June 

Statistics, time interval Feb 

Strip-line design, broadband Nov 

Synthesis, frequency Sept. 

Nov 

Synthesis, speech Jan 
Stimulus/measurement units June 
SUN operating system Mar 
Supervisor board, power supply May- 
Switching matrix, test head June 
Switching, "dry" technique |une 
Symbolic interface, software 
analyzer June 



Tasks Aug. 

Temperature control. LC column Apr 

Terminal operating system (TOSJ Aug 

Terminal, industrial Nov 

Terminals, character set differences Sept 
Terminology, board test 

connections Oct. 
Test instruction set, 4062A Test 

System June 

Test station, power supply May 
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Test system, in-circuit board Oct 

Test system, parametric June 
Testing, automatic program 

generation Oct. 

Testing. EMI. computer May 
Testing, environmental, computer May 

Testing, environmental. TIMS Sept 

Testing, open-circuit Oct. 

Testing, safeguards Oct. 

Testing, short-circuit Oct. 
Test programming, interpreter-based Oct. 
Thermal printer, personal computer Aug. 

Thinkjet Mar. 

Throughput, digital testing Oct 

Time capture measurements Dec. 

Time interval averaging June 

Time interval measurements June 

Timing diagram, logic, analysis Feb. 

TIMS Sept. 

Toleranced design May 



Touchscreen personal computer 


Aug. 


Track follower, disc drive 


Ian. 


Transactions, disc drives 


Ian 


Transient analysis 


Dec. 


Translation, software 


Sept. 


Transmission Impairment 




Measuring Set 


Sept. 


Transmitter, TIMS 


Sept. 


Tree, linguistic rules 


Sept. 


Triggers, counter, main and 




delaved 


June 


Two-counter time interval 




measurement 


June 


Twisted-pair test wiring 


Oct. 


u 




Universal counter, gated 


June 


UNIX, implementation, HP 9000 




Computers 


Mar. 


User microprogramming 


Feb. 



UWis detectors. LC Apr. 
V 

Valve, rotary. LC Apr. 

VCL. vector control language Oct. 

Vector math, network analysis Nov. 

Vibration measurements Dec. 

Virtual control panel Feb. 

Virtual front panel June 

VLFM code Jan. 
Voice circuit testing. Bell-standard Sept. 

W 

Waveform measurement library June 

Wire frame. TIMS Sept. 

Word processing software Aug. 

Workstation, engineering Mar. 

May 

z 

Zoom operation, digital filter Dec. 
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HP-71B 


Handheld Computer 


July 


HP 9000 


Model 530 System Processing 




HP Series 80 


Computers 


Jan. 




Unit 


May 






Apr. 


HP 9000 


Model 340 System rrocessing 




HP 150 


Personal Computer 


Aug. 




Unit 


May 


A600 


Computers 


Feb. 


HP 9000 Option G02 


Detached Keyboard. 




A700 


Computers 


ret). 


HP 9000 Model o20(Jomputer 


May 


A900 


Computers 


reb. 


16320A 


Pin Board 


June 


HP 1000 A-Sertcs 


Computers 


Feb, 


19800A/B 


Waveform Measurement Library 


June 


1 1 Lb -1 ... 1 

HP 1000 


Computer 


Jan. 


19801A/B 


111 f LV t 4 T i Lv — 'm H* r 

Waveform Measurement Library 


June 






Oct. 


27201A 


Speech Output Module 


|an. 


1040A 


HPLC Detection System 


Apr. 


27203A 


Speech Library 


Inn 

Jan. 


1090 


Liquid Chromatograph 


Apr. 


27205A 


Speech Library 


)<in. 


1.J4SA 


Digital Graphics Display 


Nov. 


29340s 


oingie Din oneei reeaer 


Nov. 


1965A 


Gated Universal Counter 


June 


35677A/B 


S-ParameterTest Sets 


Nov. 


1980A/B 


( )m i I lose ope Measurement 




45405A 


VisiCalc*150 


Aug. 




System 


June 


4541 OA 


Scries 100/Graphics 


Aug. 


2199C/D.2439.2139A 


A900 Computer 


Feb. 


45420A 


MemoMaker 


Aug. 


2197C/D.2437.2137A 


A700 Computer 




45422A 


Personal Card File 


Aug. 


2107AK 




Feb. 


45423A 


Financial Calculator 


Aug. 


2196C/D.2436.2136A 


A600 Computer 




46010 


Keyboard 


Aug. 


2106AK 




Feb. 


64000 


Logic Development System 


Feb. 


2G74A 


Thermal Printer 


Aug. 






June 


2932A 


General Purpose Printer 


Nov. 


6431 OA 


Software Performance Analyzer 


June 


2933A 


Factory Data Printer 


Nov. 


64600S 


Logic Timing Analyzer 


Feb. 


2934A 


Office Printer 


Nov. 


79835A 


Solvent Delivery System 


Apr. 


HP 3000 


Computer 


Jan. 


79846A 


Autoinjeclor 


Apr. 


3065 


Board Test System 


Oct. 


79847A 


Autosampler 


Apr. 


3065C 


Controller 


Oct. 


79880A 


HPLC Detection System 


Apr. 


3065H 


Test Station 


Oct. 


79881A 


Filterphotometric Detector 


Apr. 


3081 A 


Industrial Workstation Terminal 


Nov. 


82400A 


Card Reader Module 


July 


3561A 


Dynamic Signal Analyzer 


Dec. 


82401 A 


HP-1L Interface Module 


July 


3577A 


Network Analyzer 


Nov. 


82420A 


4K Memory Module 


July 


4062A 


Semiconductor Parametric Test 




82441A 


FORTH/AssemblerPac 


July 




System 


June 


824HOA 


Math Pac 


July 


4084A 


Switching Matrix Controller 


June 


82484A 


Curve FittingPac 


July 


4085A 


Switching Matrix 


June 


82967A 


Speech Synthesis Module. 




4141A 


DC Source/Monitor 


lime 




Series 80 


Jan, 


428UA 


l-M Hz CMeter/C-V Plotter 


lutie 


92188 


Character Fonts 


Nov. 


4945A 


Transmission Impairment 




97050A 


BASIC Language System. 






Measuring Set 


Sept. 




HP 9000 Model 520 


Mar. 


7933 


Disc Drive 


Jan. 


97052A 


3D Graphics. HP 9000 




7935 


Disc Drive 


Jan. 




Model 520 


May 


HP 9000 


Series 200 Computers 


Mar. 


97070A 


HP-UX, HP 9000 




HP 9000 


Serins 500 Computers 


Mar. 




Model 520, single-user 


Mar. 






May 


97079A 


HP-UX, HP 9000, Models 




HP 9000 


Model 520 Computer 


May 




530 and 540. single-user 


Mar. 
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97080A 
97089A 
97935A 



HP-UX, HP 9000 

Model 520. multiuser 
HP-UX, HP 9000 Models 

530 and 540. multiuser 
Data Pack 



Mar. 

Mar. 
)an. 



98760A 
98770A 
98780A 



Color Monitor Assembly May 

High-Performance Color Display May 
Monochromatic Monitor 

Assembly May 



PART 4: Author Index 



Abell. Stephen 


July 


Epstein. James S. 


Dec. 


Akiyama. Tomoyuki 


June 






Allan. Marlu E. 


Feb. 


Fay. Thomas R. 


Oct. 


Allen, Kevin W. 


May 


Febvre, Paul F. 


May 


Anderson, Steven R. 


May 


Felsenthal. Howell R. 


May 


Arkin. Karlie J. 


Aug. 


Firdmann. Jacques 


Nov. 






Fotland, David A. 


Feb. 


Baker, Alan |. 


Nov, 


Fritz. Gary D, 


May 


Balliew, Robert E. 


Oct. 






Balza, |ohn I. 


Mar. 


Galen. Peter M. 


Jan. 


Banno, Takuo 


Ju lie 


Geiger, Wolfgang 


Apr. 


Barbera, Joseph D. 


Aug. 


Georg, Dennis D. 


Mar. 


Batey, Robert M. 


Jan. 


Glatz. Bernd 


Apr. 


Becker. James D. 


Jan. 


Goldsmith, Kurt R. 


Sept. 


Bell, R. Frank 


Jan. 


Grodd, Laurence W. 


July 


Blasciak, Andrew J. 


June 






Blascow, Stanley M., Jr. 


July 


Hadbawnik, Detlev 


Apr. 


Bohley. Thomas K. 


June 


Hamilton. Gail E. 


I u ne 


Bounaix. Jean 


Nov. 


Hammond, Donald L. 


Mar. 


Bovven. Michael K. 


May 


Hancock. Johnnie L. 


June 


Brooks, Robert L. 


May 


Hartl, Hans-Georg 


Apr. 


Brown. Patricia A. 


Aug. 


Hartman, Douglas O. 


Feb. 


Bui. Xuan 


May 


Harwood, Vance R. 


Oct. 


Burkman. Jack L. 


May 


Hastings, Ernie 


Sept. 


Bury, Robert J. 


Mar. 


Hawk, Joseph A. 


June 




May 


Hendricks, T. Michael 


Oct. 






Hetrick. Michael V. 


Mar. 


Carlson. John R. 


Feb. 


Hill, Charles P. 


Sept. 


Carrie, Susan E. 


Aug. 


Hiller, Donald R. 


Dec. 


Carver. Brett K. 


June 


Holmberg. Randy W. 


Oct. 


Chatlerton, Craig B. 


Feb. 


Hoog, Bryan C. 


Dec. 


Christofanelli, Paul C. 


May 


Hoschele. Giinter 


Apr. 


Chu, Billy 


Feb. 


Hueftle. Elizabeth R. 


Jan. 


Clark. Elizabeth A. 


Feb. 


Hunt. Jackie 


July 


Connor, Michael L 


Mar. 






Cooley, Jack D. 


Mar. 


Ishiguro, Kenzo 


|une 


Cross, Don M. 


Feb. 






Culbertson. W. Bruce 
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Custom Digital Filters for Dynamic Signal 
Analysis 

by Charles R. Panek and Steven F. Kator 



DIGITAL FILTERING offers some powerful advan- 
tages when used with a fast Fourier transform (FFT), 
because it can perform frequency domain analysis 
in a narrow band around some arbitrary center frequency. 
To understand the need for a digital filter in dynamic signal 
analysis, let's first review the basics of FFT analysis and 
sampled systems. 

Basic FFT Analysis 

The fast Fourier transform is a special case of the discrete 
Fourier transform, a mathematical algorithm that translates 
a series of samples of a signal, evenly spaced in time, into 
a similar series of samples of the frequency spectrum of 
the signal. A finite set of these evenly spaced time samples 
is known as a time record. A time record with 1024 samples 
spaced 100 microseconds apart will be transformed into a 
frequency series or spectrum of 1024 evenly spaced sam- 
ples between dc and the 10-kHz sampling frequency f s . 
The frequency domain information for frequencies above 
iJ2 is redundant, that is, the spectrum is essentially folded 
around f s /2. This phenomenon determines the minimum 
rate at which a signal of a certain bandwidth can be sam- 
pled. Thus, provisions must be made to ensure that the 
sampled signal contains no frequency components of sig- 
nificance higher than f 5 /2 so that the FFT results do not 
fold back on themselves, an effect called aliasing. 

A very basic FFT-based signal analyzer (Fig. 1) consists 
of a sampler preceded by a low-pass antialiasing filter with 
a cutoff frequency less than f s /2, and an FFT processor that 
operates on fixed-length time records. But what if the user 
has an interest in signals at very low frequencies, or in 
signals spaced very close together in frequency? With a 
sampling rate as described above, this analyzer would only 
have about 10 Hz of frequency resolution. How can this be 
improved? One possibility is to lengthen the time record. 
A 10,000-point time record would allow 1-Hz resolution. 
However, it would take 10 times as long to capture the 
data, and roughly 13 times as long to calculate the FFT. 

Another possibility is to reduce the sampling frequency. 
This reduces the spacing of frequency components by a 
similar amount; a 1-kHz sampling frequency gives about 
1-Hz resolution with a 1024-point transform. This would 



be an ideal solution, except that the cutoff frequency of 
the low-pass input filter must then be changed to ensure 
that the spectral content of the sampled signal is entirely 
below half the new sampling frequency. To allow for many 
different bandwidths, we would have to be able to switch 
between many analog low-pass filters, or perform a compli- 
cated retuning of a single filter. 

There is another way to lower the sampling rate, how- 
ever, and that is by reducing the bandwidth of the signal 
after it is sampled and then discarding samples at regular 
intervals. This process is known as sample rate decimation, 
and the low-pass filtering is performed by a digital decima- 
tion filter. To increase the resolution by a factor of five, 
we would pass the output of our sampler through the digital 
filter, filtering the data to one-fifth the original bandwidth, 
and then discard four out of every five samples. We can 
do this because the filtered signal is now oversampled, that 
is, it contains redundant information. 

Digital filtering offers some other special advantages. The 
performance of a digital filter is a function of the coefficients 
used and the length of the data word. Once a digital filter 
design is completed, the bandwidth, flatness, noise, phase 
response, and other parameters of the filter are exactly 
repeatable from circuit to circuit. The adjustments and toler- 
ance variances characteristic of analog filters are eliminated. 

Zoom Operation 

By cascading decimating digital filters, we can reduce 
the bandwidth and increase the frequency resolution by 
larger amounts. We still have an important restriction, how- 
ever: the maximum frequency of our analysis continually 
drops lower. Suppose we wish to look at the modulation 
sidebands of a 4-kHz carrier frequency. If the modulation 
is very low in frequency, say 1 Hz, we will be unable to 
resolve the sidebands with the simple system of Fig. 1. 
Simply decimating the converter output by a factor of two 
will filter out the carrier, sidebands and all, assuming a 
10-kHz sampling frequency. 

The nature of the FFT limits us to frequency analysis of 
a time series whose frequency content extends from dc to 
half the sampling frequency. We can overcome this limita- 
tion by borrowing a technique from the designers of radio 
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receivers and swept-frequency spectrum analyzers. That 
is. demodulation or analysis is performed in a fixed inter- 
mediate frequency (IF) band, and any frequency band of 
interest is mixed with a local oscillator (LO) signal to trans- 
late it to the IF band. 

Here, we can consider our IF band to extend from dc up 
to half the sampling frequency, adjustable in bandwidth 
by the decimation rate of the digital filter. Thus, we can 
mix the 4-kHz carrier frequency with a 4-kHz local oscil- 
lator frequency and analyze the dc and low-frequency prod- 
ucts using a decimated sample rate that will give us the 
desired resolution. In addition, we can perform quadrature 
mixing, that is. multiply the input by sine and cosine waves 
at the LO frequency. By processing the resulting complex 
numbers in a complex FFT. frequency data both above and 
below the LO frequency can be obtained without generating 
unwanted images. We can once again use the advantages 
of digital signal processing to do this mixing in a stable 
and accurate way. The output of the analog-to-digital con- 
verter (ADC) is fed into a pair of digital multipliers, along 
with sine and cosine wave sequences of the desired center 
frequency sampled at the same rate. The resulting digital 
products are fed into a pair of identical digital filters for 
sample rate reduction. This is referred to as "zoom" oper- 
ation, for it allows the user to expand the view of a given 
portion of the input signal spectrum. 

Architecture and Algorithm 

We have arrived at the basic digital signal processing 
block diagram of the HP 3561 A Dynamic Signal Analyzer 
(Fig. 2). A 100-kHz analog antialiasing filter precedes a 
sample-and-hold circuit and ADC running at a sample rate 
of 256 kHz. This sample rate gives nicely spaced frequency 
points when used with 2 N -sized time records (in this case. 
1024 samples), The digital LO generates cosine and sine 
sequences, sampled at the ADC sample rate, which are fed 
into the mixer multipliers of the real and imaginary com- 
ponent digital filters, respectively, along with the ADC out- 
put. In baseband mode, where no frequency translation is 

Signal 
Input 

V 



desired, the cosine output of the LO is a constant, and only 
the real component filter is used. The 20-bit output of the 
filter chips is rounded to 16 bits by the control 1C and fed 
into RAM, where it is accessed by the FFT processor. 

It was determined on the basis of noise performance and 
data scaling that the digital filters should use 20-bit num- 
bers to represent internal data. The A-to-D sample rate is 
specified at 256.000 samples per second. Hence, to allow 
serial processing of a decimated sample between each new 
ADC sample, a clock rate of 2 (samples) x 20 (bits) x 256 
kHz. or 10.24 MHz. would be required. 

Because of power limitations and other concerns, an al- 
ternative to strictly serial processing is used. Bits of the 
data word are processed in pairs of adjacent bits. Thus, the 
two least-significant bits (LSBs) are processed on the first 
clock cycle, the next two bits are processed on the second 
cycle, and so forth. This lowers the required clock rate to 
5.12 MHz. 

Mixer Multiplier 

The multiplier used in the HP 3561A's signal processing 
circuitry uses a serial paired-bit shift-and-add scheme that 
operates on two's-complement words. To understand its 
operation, consider a basic serial/parallel unsigned multi- 
plier as shown in Fig. 3. 

The circuit is initialized by clearing all the D flip-flops 
except the one flip-flop at the rounding point, which is 
set. and then loading the multiplicand into the parallel 
latch. As each multiplier data bit is brought in (least-signif- 
icant bit first), a new partial product is formed at the AND- 
gate outputs. The partial products are accumulated by the 
carry-save adders, and the result is presented serially at 
the output, with one new product bit for every new multi- 
plier bit. The circuit is clocked until all the product bits are 
shifted out, and the input data is padded with zeros as 
necessary. 

We may enhance this design to accommodate two's-com- 
plement numbers by adding a serial negator ahead of the 
most-significant AND gate in the serial input path. If the 
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multiplicand is negative, this has the same effect as sub- 
tracting the last partial product. 

This design was adapted to work with bit pairs by split- 
ling the serial inpul data stream into odd and even bits, 
and generating three two-bit-wide data streams represent- 
ing the data multiplied by 1. 2. and 3. The parallel data 
(multiplicand) is also grouped into pairs. The AND gates of 
the basic circuit are replaced with pairs of 4:1 multiplexers, 
with the multiplicand bit-pairs driving the select lines and 
the XI, x2.and x 3 multiplierbit-pairsdrivingtheinputs. 

The final design is shown in Fig. 4. The ADC and LO 
data words are brought in serially and converted to bit 
pairs. The LO word enters a two-bit-wide, serial-in. paral- 
lel-out shift register. The ADC word passes through a hard 
limiter that prevents overflowed ADC data from entering 
the filter, and a 0.75 scaler, which reduces the amplitude 
of the data so that overshoot in the filter will not result in 
overflows. This multiplier outputs the 20 most-significant 
bits (MSBs) of the product of a 16-bit LO word and a 13-bit 
ADC data word in ten clock cycles, which suits the require- 
ments of the digital filter perfectly. 

Cascaded, Pipelined, Second-Order Sections 

There are two filters on the HP 3561A's custom digital 
filter chip. One. with a bandwidth of approximately 20% 
of the input sample rate, is used before reducing f s by a 
factor of 2. The other filter, with a 7.8% bandwidth, is used 
before decimation by a factor of 5. Both of these filters are 
implemented as a cascade of second-order sections. Both 
filters also share the same digital circuitry. Filter selection 
is accomplished by an additional address bit to the on-chip 
ROMs that contain the filter coefficients. Each filter is an 
eighth-order elliptical design, with coefficients chosen to 
meet design criteria while minimizing the number of bits 
to represent them. Having fewer bits in the coefficients 
reduces the amount of digital hardware needed to imple- 
ment the filters. 

Pipelining of the circuitry is used extensively between 
the second-order sections of the filter. That is, data is de- 
layed by one sample between each filter section. Thus, the 
entire filter processes several samples at one time, which 
increases throughput without increasing clock frequency. 

There are ten channels of data storage associated with 
each filter section and the output of the filter. Using this 
memory, data can be recycled through the filter up to ten 
times after the sample rate has been reduced each time. By 
processing decimated data between samples from the ADC, 
the filter sample rate need only be twice that of the A-to-D 
converter. 1 

Peled-Liu Second-Order Sections 

Each second-order filter section must implement a sim- 
ple difference equation: 

Yn = 3oX n + a^., + a 2 x n 2 - b,y n _, - b 2 y n _ 2 (1) 

where x n is the nth input word, y M is the nth output word, 
and a, and b| are the filter coefficients. 

To implement this equation with a minimum of 
hardware, a technique variously known as the Peled-Liu 2 
or Princeton multiplier is used. A brief explanation of this 



technique and its adaptation to this design follows. If we 
assume that data is represented by two's-complement bi- 
nary numbers of B bits bounded by 1 and -1. we can 
represent x n of Equation 1 by: 

B- 1 

*n=-xS+2 xi2~' x» = 0 or 1 (2) 



where x, 1 , is the jth bit of the nth input word. A similar 
expression represents y„ of Equation 1. Using these expres- 
sions, the difference equation can be rewritten as: 




We can define a function tli with five binary arguments as: 
iMx', x-. x 3 , x\ x 5 )= a 0 x l + a 1 x 2 + a 2 x 3 -b l x 4 -b z x s (4) 

where x k = 0 or 1 
Then Equation 2 can be rewritten as: 

H-l 

v„ = ^ 2~> iHxl.xf, -LX^.y/.^.y, 1 ,-;,) 
i=i 

-■Mx' n ',x!! „xJUy°_„y°_ 2 ) (5) 
The function i//can have only 32 possible values, depend- 




Fig. 5. 8as/c Peled-Uu second-order tiller Blocks A, B, and 
D are serial-in, serial-out shitt registers. Block C is a parallel-in, 
serial-out shitt register 
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ing on its five binary arguments. Since the values of a> and 
bj are fixed in each filter section, this suggests implement- 
ing this function in ROM. The entire second-order section 
is shown in Fig. 5. 

Shift registers A through D are all B bits long. The nth 
input word x„ is fed, least-significant-bit first, into the ROM 
and shift register A. The ROM is simultaneously addressed 
by the corresponding bits of the x„_,, x n _ 2 , y„-i and y„_ 2 
words. The ip value generated on each clock cycle is ac- 
cumulated with the result of the previous operation, shifted 
right one bit. On the first cycle, the accumulator register 
is preloaded with a rounding value, which is added to 
assure minimum bias when the result is truncated back to 
B bits. On the final cycle of each filter operation, the t/i 
value corresponding to the sign bits of the data word is 
subtracted from the accumulator, and the accumulator 
value is parallel-loaded into shift register C. 



Paired-Bit Implementation 

To maintain a lower clock rate, as mentioned earlier, the 
design of Fig. 5 was adapted to process bit pairs, using two 
Princeton multipliers per second-order section. Equation 
5 can be rewritten: 

B 1 

H- 2 

+ 2 2 "V (xA.xA-,.xA_ 2 .yi_ 1 .y» 2 ) 

1=2.4.... 

-#&j^3<§ 2 .yn- 1 .yn- 2 ) (6) 

where B is an even number. 
This leads to a second-order section design as shown in 
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Fig. 6. Basic circuit of Fig. 5 
adapted to process bit pairs using 
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Fig. 6. Memory is organized as banks of 10-bit shift registers. 
One bank stores the even bits of the 20-bit data words, and 
the other stores the odd bits. Each bank consists of 10 
channels of memory, one for each pass the data can make 
through the filter. Thus, a second-order section can filter 
one sample from the output of the mixer multiplier, and 
then store away its state and operate on a sample from a 
data sequence that has already passed through the filter 
and had its sample rate reduced. 

Each Princeton multiplier generates a partial product. 
20 bits long, after rounding. These numbers are then added, 
two bits at a time, in the carry/save adder. The resulting 
sum is then split into its odd and even bits, and fed into 
the appropriate channel of y n _, memory. 

The y n _ , memory for the first filter section serves double 
duty as the x n _, memory for the second section. Pipeline 
registers are used for delaying and resynchronizing the data 
between sections. 

The simplified block diagram for the complete filter is 
given in Fig. 7. Note that actually five second-order sections 
are used to implement an eighth-order filter. The first sec- 
tion implements a transfer function consisting of a scaling 
constant and two poles, and the last section implements a 
transfer function with only a pair of zeros. This implemen- 
tation was chosen to minimize circuit area on the IC. 

The last section of the filter is followed by an overflow 
detect and hard-limit circuit. Although the data is scaled 
down at the input so that internal overflows cannot occur, 
to preserve unity gain, the data is scaled up in the last 
section. If an overflow does occur there, the output data is 
set to the appropriate ± full-scale value. 

Once data emerges from the hard limiter, it is stored in 
the 10-channel-by-20-bit feedback memory, where it is re- 
circulated until new data is written into the same channels. 
Thus, sample-rate reduction of filtered data is accom- 
plished by merely writing into feedback memory more often 
than reading from it. 

Limit Cycles 

fust as there are challenges that face the designer of an 
analog filter who wishes to maximize noise performance 
and dynamic range, the digital filter designer must cope 
with problems unique to digital signal processing. These 
include noise-like signals and dc offsets generated within 
the filter because of errors in rounding, and low-level self- 
oscillations in the absence of a signal. 

The phenomenon of limit cycle oscillations occurs in 
many digital signal processing systems where there is a 
rounding operation within a feedback loop. As a simple 
example, consider the following difference equation: 



n>0 



y„ = x n + 0.9 y,,., (7) 

If we set y n _ , = 0 and x„ = 1 when n = 0. and x n = 0 when 
>0. then the following sequence results: 



y 0 = i 

y, = 0.9 

y, = 0.81 

y 3 = 0.729 

y 4 = 0.6561 



As can be seen, the sequence decays exponentially to- 
ward zero. Now let's limit the precision of the arithmetic, 
and round each result to the nearest tenth. Then: 



y<> 


= i 




yt 


= 0.9- 


*0.9 


y 2 


= 0.81 


— 0,8 


y 3 


= 0.72 


~*0.7 


y* 


= 0.63 


-►0.6 


ys 


= 0.54 


~*0.5 


ya 


= 0.45 


~-*0.5 


y> 


= 0.45 


-*0.5 



This time the sequence "hangs up" at 0.5; even with no 
input there is an output. A dc output like this is known as 
a deadband. When there are two feedback terms in our 
difference equation as in Equation 1 (thus, implementing 
a pole pair), the output can decay into an ac oscillation 
known as a zero-input limit cycle. Bounds for the 
amplitudes of these limit cycles have been developed. 3 In 
general, they relate to the Q of the poles of the transfer 
function. As the Q of the poles increases, limit cycles of 
higher amplitude become more likely. Since sharp-cutoff 
filters tend to have high-Q poles, a digital filter designer 
can expect to have problems with limit cycles when trying 
to maximize performance. 

When incorporating a digital filter in an FFT analyzer, 
the presence of limit-cycle oscillations is especially unde- 
sirable. Because these oscillations tend to be roughly 
sinusoidal near the resonant frequency of the pole pair, 
they will appear as a spurious signal in the frequency dis- 
play of the analyzer when the input is removed or held at 
very low levels. To add to the confusion, the oscillations 
can have different values or go away altogether, depending 
on the trajectory the Filter input takes as it decays. 

Limit cycles are not readily predictable in complicated 
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Fig. 7. Simplilied block diagram of 
complete digital filter. 
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fillers, For a basic second-order filler, it is a simple matter 
to find all the zero-input limit cycles by repeatedly allowing 
the filter output to decay (or settle into oscillation) after 
presetting the filter's registers to every possible combina- 
tion of small values (one pair at a time). However, the 
results of such an experiment do not guarantee the limit 
cycle performance when the filter input is other than zero, 
or when the filter is cascaded with other second-order sec- 
tions. 

Deadbands and limit-cycle oscillations can be explained 
qualitatively as the rounding operation dominating the fil- 
ter coefficients in determining the decay of the filter's out- 
put. A design tactic that has often proven successful is to 
inject a low-level dither signal at the rounding point. In 
the example shown above, subtracting a value of 0.05 from 
the number before rounding will allow the filter output lo 
drop below 0.5. If we add 0.05 as often as we subtract it, 
the rounding operation will still yield a minimum dc offset. 

This solution has its disadvantages, however. If a simple 
periodic dither signal is used, it too will appear as a spuri- 
ous signal on the output of the FFT stage, albeit of predict- 
able amplitude and frequency. As a matter of practice, such 
a dither signal often must be of an amplitude equal to or 
greater than the amplitude of the oscillations it must break up. 

In the HP 3561 A Dynamic Signal Analyzer, limit cycles 
are eliminated by using a dither signal derived from a long 
pseudorandom noise generator. This proves quite effective 
in eliminating deadbands and limit cycles altogether. The 
energy from this noise generator is not concentrated at one 
frequency, but is spread evenly over the entire spectrum 
when its sequence length is sufficiently long. Thus, instead 
of a spurious component at some frequency on the FFT 
display, the dither appears as an increase in the noise floor. 
In fact, the peak value of the pseudorandom noise dither 
at each rounding point in the filter is less than one LSB of 
the 20-bit data word. Therefore, not only are there no limit 
cycles or deadbands, but the dither generator is not visible 
at the output under zero input conditions. 

Offsets 

Offsets or errors as a result of rounding can result in 
significant measurement inaccuracies. Offsets in the filter 
output are most problematic during zoom measurements, 
because input signals at the same frequency as the local 
oscillator are dovvnconverted to dc and are displayed at 
midscreen. A dc error in the filter output thus results in a 
spurious component at the center frequency on the display. 

One would expect that since the filter data is 20 bits 
long, any reasonable rounding scheme should result in 
acceptably low distortion. However, several effects exacer- 
bate the problem. 

Since the HP 3561A filter is a cascade of five sections, 
and since the decimation algorithm entails further cascad- 
ing the filter with itself up to ten times, any consistent 
offset in the section arithmetic will accumulate. Also, since 
rounding occurs after each multiply within the second- 
order sections, and because of the quantization of the coef- 
ficients, the results of certain coefficient multiplications 
have nonuniform statistical distributions. Hence, biases 
occur. Finally, the dc gain from some internal nodes to the 
output of the filter can be fairly high, particularly with the 



higher-Q, divide-by-5 coefficients. 

Curing DC Offsets 

As a first example, consider the problem of rounding 20 
bits to 16 bits. If we simply truncate the unwanted bits, 
the magnitude of the error introduced will be distributed 
between 0 and - 15 LSBs of the 20-bit word. If we assume 
that the value of the truncated bits is uniformly distributed, 
the error will have a mean of -7.5 LSBs (see Fig. 8a). 

The traditional solution is to add a value of 8 LSBs \Vi 
LSB of the 16-bit result), which yields a mean error of VS 
LSB of the 20-bit word. This is quite good, but not good 
enough for our application. 

Also consider the case of nonuniformly distributed val- 
ues that are truncated. For example, suppose the first bit 
below the truncation point has a 75% probability of being 
zero. The truncation error will be distributed as in Fig. 8b, 
with a mean of - 5.5 LSBs. If we do the standard procedure 
of adding 8 LSBs. an error of 2.5 LSBs still exists. 

The solution employed in the HP 3561A digital filter 
consists of adding one of Iwo predetermined values before 
truncation. The two values are calculated so that their mean 
added to the mean truncation error equals zero. The choice 
of which value to add is controlled by a pseudorandom 
register bit with a 50% duty cycle. Thus, on the average, 
the rounding error will be zero. This technique has two 
major advantages. It is flexible enough to deal with highly 
nonuniform distributions, and it allows the designer lo 
adjust the variance in the rounding (noise power) to suit 
such needs as breaking up limit cycles. 

To see how this so-called "dithered rounding" is applied, 
return lo the first example above (Fig. 8a). If we alternately 
add zero and 15 LSBs to the value before truncation, there 
will be an average error of zero in the operation. Likewise 
in the second example (Fig. 8b), we may alternately add 
zero and 11 LSBs before truncating, again achieving zero 
mean error. Clearly, this has profound effects on the round- 
ing noise, but the average noise power will generally be 
less than four times that of simple truncation. The use of 
random alternation is important, because the alternative 
use of periodically determined rounding is equivalent to 
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Fig. 8. Example of probability mass function oi truncation error 
for (a) uniformly distributed values and (b) nonuniformly distrib- 
uted values of LSBs when truncating a20-bit word to 7 6 o//s. 
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modulating the signal, which is not normally desirable in 
a digital filter. 

Mask Design Process 

The digital decimation filters used in the HP 3561A are 
implemented as a set of three custom NMOS integrated 
circuits, fabricated at HP's Loveland Technology Center. 
Two identical filter chips, containing 30.000 transistors 
each, perform simultaneous filtering of the real and imag- 
inary data. A control chip containing about 10.000 transis- 
tors provides correct sequencing of control signals for the 
filter ICs and provides an interface between them and the 
main instrument processor and the memory used by the 
FFT processor. 

Designing integrated circuits of this complexity requires 
a methodology that is logical and well thought out, and 
that leaves no room for errors. Once the filter algorithms 
and topology are determined, and the digital design com- 
pleted, the process of translating schematics into silicon 
begins. This process starts with preliminary layouts of crit- 
ical portions of the circuitry, both to help prove feasibility 
of the design, and to gain insight into the area requirements 
of the circuit. A large portion of the circuitry in the digital 
filter consists of shift register memory. A simple shift regis- 
ter cell layout was optimized for size, and, along with a 
few other cells, this critical cell was used to develop an 
architectural plan showing the sizes and locations of the 
major blocks of circuitry and the power and signal routing 
between these blocks. 

The artwork for the chips is built in a hierarchical fash- 
ion. The chip is divided into large functional blocks, (mul- 
tiplier, filter arithmetic, memory, et cetera). These blocks 
are rectangular and their horizontal and vertical dimen- 
sions match. The large blocks are themselves composed of 
smaller, regularly shaped blocks of medium complexity 
(control circuitry, parallel arithmetic, input latches, etcet- 
era). This hierarchy continues down to the most basic 
blocks of circuitry, called leaf cells (e.g.. shift registers and 




(a) 



flip-flops). 

A keystone of this design method is the use of a standard 
metal pattern. The metal layer on the IC is used to distribute 
power and route signals over long distances. Power distri- 
bution on the chip is done using a metal pattern than looks 
like a pair of interdigitated combs (Fig. 9a). One set of 
fingers is connected to the 5V supply, the other to ground. 
Thus, the ground and 5V lines are spaced at regular inter- 
vals, alternating across the chip. Another element of the 
standard metal pattern is that room for a few metal signal 
lines is left between the power lines. 

All leaf cells are built using the standard metal pattern 
(Fig. 9b). Every cell must be the same height, and includes 
half of a 5V line at the top and bottom of the cell, and a 
ground line across the middle of the cell. Between the 
power lines is room for up to four metal signal bus lines. 
Circuitry is laid out underneath these metal lines, using 
either the metal bus lines or portions of the underlying 
polysilicon layer for interconnect. A half-height cell can 
also be used for smaller circuit blocks. Thus, chip area that 
might otherwise be used just for busing signals from one 
block to another can also contain circuitry. Care must be 
exercised in layout, of course, to ensure that noise-sensitive 
nodes are not placed under clock or data lines. 

Production Considerations 

It is not sufficient that an IC design be functional in an 
instrument under nominal conditions. Design goals for the 
chip must be developed early that will permit reliable op- 
eration given variances in ambient temperature, supply 
voltages, fabrication processing, and the performance of 
external circuitry. 

Testing. Given these criteria, production testing of these 
integrated circuits is a two-part process. Each IC must be 
thoroughly checked for correct functional operation of all 
of its many thousands of transistors, then this operation 
must bi> confirmed under test conditions that attempt to sim- 
ulate worst-case conditions of operation in the instrument. 



VDD 




^^^^^^^^^^ Fig. 9. (a) Simplified layout ol top 
metal pattern tor digital tiller chip. 
^^L^S (b) Standard metal pattern used in 

designing cells lor the digital filter 
(b) chip (see shaded area in (a)) 
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Complete functional testing of large, complex chips pre- 
sents a special challenge, especially with a circuit like the 
digital filter, whose data inputs are all serial. Determining 
a sequence of ones and zeros on these inputs that is guaran- 
teed to excite all of the many thousands of nodes on the 
chip is a virtually impossihle task, given the length of time 
a complete simulation of the entire chip would take. 

Because most of the internal data path of the digital filter 
is also serial, direct accessibility of internal nodes is not 
the problem it might otherwise be. If any node in the serial 
data path is stuck at a one or zero, this condition will 
propagate an error to the output. Other, nonserial portions 
of the circuit, such as ROMs, parallel adders, and control 
circuits, are tested by operating the chip in each of its 
modes, and exciting the filter with signals of large dynamic 
range for as long a test as is feasible. In addition, a special 
test mode allows use of any of the 10 channels of data 
memory without passing data through the previous chan- 
nels first, thus speeding up the memory test and helping 
to isolate any problems to a specific channel. 
Test Program Generation. Rather than manually generate 
a set of stimulus and response vectors to test the chips and 
then code these vectors into test system programming lan- 
guage one clock cycle at a time, a software system was 
developed that automatically generates test system lan- 
guage code when given a stimulus sequence. This system 
generates response vectors using the same software 
simulator that was used to model the chips during the 
design process. Other software then converts these vectors 
into serial data streams used by the chips and generates 
digital test code for the wafer and packaged parts test sys- 
tem used in production. 

Packaging, Cooling 

Dense, high-speed chips like the digital filter and control 
chips require special consideration for packaging and cool- 
ing. A 64-pin. grid-array package was chosen for both parts. 
This choice not only provides the needed number of pins 
for the control chip, but offers a low thermal resistance 
between the substrate and the top surface of the package. 

A cast aluminum heat sink that fits over the chips on the 
printed circuit board in the HP 3561 A Analyzer promotes 



cool, reliable operation even under conditions of elevated 
ambient temperature. The heat sink mounts to the board 
using a spring loading scheme that keeps the heat sink's 
mating surface flush and in firm contact with the surface of 
the package without unduly stressing the brittle ceramic 
package. 
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